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ABSTRACT 



In this paper we present a clustering analysis of QSOs using over 20000 ob- 
jects from the final catalogue of the 2dF QSO Redshift Survey (2QZ), measuring 

■ the redshift-space two-point correlation function, £(s). When averaged over the red- 
f^ i shift range 0.3 < z < 2.2 we find that £(s) is flat on small scales, steepening on scales 

above - 25 h' 1 Mpc. In a WMAP/2dF cosmology (fi m = 0.27, Q A = 0.73) we find 

■ a best fit power law with so = 5.48ig'4g hr 1 Mpc and 7 = 1.20 ± 0.10 on scales 
^ c| s = 1 to 25 hr 1 Mpc. We demonstrate that non-linear redshift-space distortions have 

O ^ a significant effect on the QSO £(s) at scales less than ~ 10 h^ 1 Mpc. A cold dark 

matter model assuming WMAP/2dF cosmological parameters is a good description of 
the QSO £(s) after accounting for non-linear clustering and redshift-space distortions, 
and allowing for a linear bias at the mean redshift of 6q(z = 1.35) = 2.02 ± 0.07. 
[ We subdivide the 2QZ into 10 redshift intervals with effective redshifts from z = 

0.53 to 2.48. We find a significant increase in clustering amplitude at high redshift 
in the WMAP/2dF cosmology. The QSO clustering amplitude increases with redshift 
, such that the integrated correlation function, £(s), within 20 h^ 1 Mpc is £(20, z = 

0.53) = 0.26 ± 0.08 and £(20, z = 2.48) = 0.70 ± 0.17. We derive the QSO bias and 
find it to be a strong function of redshift with 6q(z = 0.53) = 1.13 ± 0.18 and 6q(z = 
2.48) = 4.24 ± 0.53. We use these bias values to derive the mean dark matter halo 
(DMH) mass occupied by the QSOs. At all redshifts 2QZ QSOs inhabit approximately 
the same mass DMHs with Mdh = (3.0 ± 1.6) x 10 12 h^ 1 M Q , which is close to the 
characteristic mass in the Press-Schechter mass function, M*, at z — 0. These results 
imply that Lq QSOs at z ~ should be largely unbiased. If the relation between black 
hole (BH) mass and Mdh or host velocity dispersion does not evolve, then we find 
that the accretion efficiency (L/L^dd) for Lq QSOs is approximately constant with 
redshift. Thus the fading of the QSO population from z ~ 2 to ~ appears to be due to 
less massive BHs being active at low redshift. We apply different methods to estimate, 
tQ, the active lifetime of QSOs and constrain £q to be in the range 4 x 10 6 — 6 x 10 8 
years at z ~ 2. 

We test for any luminosity dependence of QSO clustering by measuring £(s) as a 
function of apparent magnitude (equivalent to luminosity relative to £q). However, 
we find no significant evidence of luminosity dependent clustering from this data set. 

Key words: galaxies: clustering - quasars: general - cosmology: observations - large- 
scale structure of Universe. 
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1 INTRODUCTION 

The question of how activity is triggered in the nucleus of 
galaxies is vital to answer if we wish to have a full under- 
standing of the galaxy formation process. It appears that 
a large fraction of galaxies may have contained an active 
galactic nuclei (AGN) at some point in their history. When 
local galaxies are surveyed (including our own Milky Way) 
most show evidence of a super-massive black hole (BH) (e.g. 
Kormendy & Richstone 1995). The BHs tend to be found 
in dynamically hot systems (i.e. spheroids - elliptical galax- 
ies or bulges), and the mass of the BHs is well correlated 
with the mass of the spheroid. The tightest correlation is 
found between BH mass, Mbh, and spheroid velocity dis- 
persion, a* (Gebhardt et al.2000; Ferrarese & Merritt 2000). 
At higher redshift it is not clear that this correlation holds, 
or indeed in general, how high redshift BHs relate to their 
host galaxies. However Shields et al. (2003) do suggest that 
the same Mbh — a seems to be appropriate at high redshift. 

It is the powerful evolution in luminosity of the AGN 
population which allows them to be readily observed to high 
redshift. Understanding this evolution goes hand-in-hand 
with our understanding of the relation between AGN and 
galaxies. Croom et al. (2004a) (which we will henceforth 
call Paper XII) find that optically selected QSOs are well 
described by so called 'pure luminosity evolution' (PLE) 
with an exponential increase in the typical luminosity Lq 
(e-folding time of ~ 2 Gyr) up to z ~ 2. Work at higher red- 
shift (e.g. Fan et al. 2001) find that at z ~ 4 — 6 the number 
density of QSOs is much lower than at z ~ 2. The X-ray 
luminosity function (LF) appears to give a more complex 
picture (Ueda et al. 2003) but still shows the general trend 
of luminous AGN being more active, peaking at z ~ 2 — 3. 

The question is then, how do we gain further informa- 
tion about the physical processes of QSO formation at high 
redshift? One approach is to attempt to directly image QSO 
host galaxies at high resolution (Kukula et al. 2001; Croom 
et al. 2004b). These analyses seem to show that high red- 
shift QSO hosts (at least for radio quiet sources) are no 
brighter than low redshift hosts, after accounting for only 
passive evolution of the stellar populations in the galaxies. 
QSO clustering measurements gives us an important second 
angle to study the hosts of QSOs, as the clustering ampli- 
tude can be considered as a surrogate for host mass or dark 
matter halo (DMH) mass, Mdh- With large samples such as 
the 2dF QSO Redshift Survey (2QZ; Paper XII) it is pos- 
sible to determine these host properties over a wide range 
in redshift. With an estimate of the host mass of these high 
redshift QSOs we can hope to determine whether the host 
mass vs. BH mass correlation at low redshift continues to 
high redshift. We can also attempt to predict the masses of 
the descendents of high redshift QSOs, and locate them in 
the local universe. 

A number of authors (e.g. Martini & Weinberg 2001; 
Haiman & Hui 2001; Kauffmann & Haehnelt 2002) have 
constructed models for QSO evolution including clustering, 
and these need to be tested against accurate measurements. 
One parameter that can be derived from these models is a 
mean QSO lifetime, although the exact interpretation of this 
is rather model dependent. 

As well as being used for the study of QSO forma- 
tion/evolution, QSOs are also powerful probes of large-scale 



structure in their own right. The large volumes probed 
(~ 6 x 10 9 h~ 3 Mpc 3 for the 2QZ in a universe with fi m = 0.3 
and Qa = 0.7) and high redshift sampled makes observa- 
tions quite complementary with lower redshift galaxy obser- 
vations and higher redshift CMB observations. A number of 
authors have attempted to detect high redshift QSO cluster- 
ing (Osmer 1981; Shaver 1984; Shanks et al. 1987; Iovino & 
Shaver 1988; Andreani & Cristiani 1992; Mo & Fang 1993; 
Shanks & Boyle 1994; Croom & Shanks 1996; La Franca 
et al. 1998) and made some preliminary measurements of 
clustering evolution, but these have all been based on small 
samples of QSOs (typically a few hundred objects). At low 
redshift, there have also been a number of recent analysis. 
Grazian et al. (2004) find so = 8.6 ±2 Mpc for a sample 
of bright, B < 15, low redshift, z < 0.3, QSOs. Miller et 
al. (2004) show that the AGN fraction in the SDSS galaxy 
survey is not dependent on environment, while Croom et 
al. (2004c) and Wake et al. (2004) show that low redshift, 
low luminosity AGN are clustered identically to non-active 
galaxies. The 2QZ provided the first large, deep sample with 
which to perform detailed clustering analysis at high red- 
shift. Outram et al. (2003), Outram et al. (2004), Miller et 
al. (2004) and others have used the 2QZ to test cosmologi- 
cal models. The two-point correlation function (the subject 
of this paper) has been discussed by Croom et al. (2001a) 
for the preliminary, 10k, data release of the 2QZ (Croom 
et al. 2001b). They found that the clustering of high red- 
shift (z ~ 1.5) QSOs to be very similar to the clustering of 
typical galaxies at low redshift. They also found that the am- 
plitude of clustering was approximately constant, or slightly 
increasing, with redshift. 

For comparison to the high redshift QSO clustering re- 
sults, there are now some measurements of galaxy cluster- 
ing over similar redshift intervals. These suggest moderately 
high clustering amplitudes, generally not inconsistent with 
that measured for QSOs. E.g. Deep wide-field (~ few de- 
grees) imaging surveys used to measure the angular correla- 
tion function of galaxies also suggest high clustering ampli- 
tudes (Postman et al. 1998). However, various differences are 
found, depending on the magnitude limits and photomet- 
ric bands used to define the samples. This is not surprising 
given that there is clearly evidence that galaxy clustering 
is a function of luminosity (Norberg et al. 2001). This may 
also be the case for QSOs, although there has been no sig- 
nificant evidence for this to date (Croom et al. 2002). At 
z ~ 3, galaxy surveys using the drop-out technique (e.g. 
Steidel et al. 1998) have found that L ~ L* galaxies also 
cluster similarly to local galaxies on scales <J10 h" 1 Mpc, 
with ro ~ 4 - 6 hT x Mpc for a cosmology with f2 m = 0.3 
and Q A = 0.7 (Adelberger et al. 1998; Foucaud et al. 2003; 
Adelberger et al. 2003). 

In this paper we use the final data relase of the 2QZ 
(Paper XII) to measure the QSO two-point correlation func- 
tion over a wide range in redshift. The 2QZ is currently the 
best sample on which to perform this analysis, being by far 
the largest QSO sample with a high surface density (~ 35 
deg -2 ). We focus in this paper on the redshift-space corre- 
lation function £(s) and attempt to account for the effects 
of any z-space distortions. The real-space correlation func- 
tion will be addressed in a further paper (da Angela et al. 
in preparation), and the cross-correlation of QSOs in differ- 
ent luminosity intervals will be discussed by Loaring et al. 
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(in preparation). In Section 2 we introduce the 2QZ sample 
and the techniques used in our analysis. In Section 3 we use 
mock QSO catalogues (Hoyle 2000) constructed from the 
large simulations to test the reliability of our corrections 
for variations in completeness in the 2QZ. The redshift av- 
eraged, redshift dependent and luminosity dependent 2QZ 
£(s) measurements are presented in Sections 4, 5 and 6 re- 
spectively. We finally discuss our conclusions in Section 7. 



2 DATA AND TECHNIQUES 

2.1 The 2dF QSO Redshift Survey 

There is a full description of the 2QZ in Paper XII. Briefly, 
the survey covers two 75° x 5° strips, one passing across 
the South Galactic Cap centred on <5 = —30° (the SGP 
strip) and the other across the North Galactic Cap centred 
on 5 = 0° (the NGP or equatorial strip). The SGP strip 
extends from a = 21 h 40 to a = 3 h 15 and the equatorial 
strip from a = 9 h 50 to a = 14 h 50 (BI950). The total survey 
area is 721.6 deg 2 , when allowance is made for regions of sky 
excised around bright stars. 

2dF spectroscopic observations were carried out on 
colour selected targets in the magnitude range 18.25 < bj < 
20.85. This resulted in the discovery of 23338 QSOs at red- 
shifts less than z ~ 3. In this paper we restrict our analysis 
to QSOs with quality 1 identifications (see Paper XII), that 
is 22655 QSOs. The distribution of QSOs in the final sample 
is shown in Fig. 1. 



2.2 Correlation function estimates 

As the QSO correlation function, £(s), probes high redshifts 
and large scales, the measured values are highly dependent 
on the assumed cosmology. In determining the comoving sep- 
aration of pairs of QSOs we choose to calculate £(s) for two 
representative cosmological models. The first uses the best 
fit cosmological parameters derived from WMAP, 2dFGRS 
and other data (Spergel et al. 2003; Percival et al. 2002) with 
(fi m ,n A ) = (0.27, 0.73), which we will call the WMAP/2dF 
cosmology. The second model assumed an Einstein-de Sit- 
ter cosmology with (O m ,f2A) = (1.0,0.0), which we denote 
as the EdS cosmology. We will quote distances in terms of 
h~ Mpc, where h is the dimensionless Hubble constant such 
that H = WOh km s _1 Mpc" 1 . 

We have used the minimum variance estimator sug- 
gested by Landy & Szalay (1993) to calculate £(s), where 
s is the redshift-space (or 2-space) separation of two QSOs 
(as opposed to r, the real-space separation). This estimator 
is 



A£( S ) 



! + €(«) 
y/QQ(s)' 



(2) 



QQ(s) - 2QR{s) + RR(s) 
RR(s) ' 



(1) 



where QQ, QR and RR are the number of QSO-QSO, QSO- 
random and random-random pairs counted at separation s± 
As/2. QR and RR are normalized to the total number of 
QSOs. The density of random points used was 50 times the 
density of QSOs. 

We calculate the errors on £(s) using the Poisson esti- 
mate of 



At small scales, <50 h^ 1 Mpc, this estimate is accurate be- 
cause each QSO pair is independent (i.e. the QSOs are not 
generally part of another pair at scales smaller than this). 
On larger scales the QSO pairs become more correlated and 
we use the approximation that A£(s) = [1 + \/Nq, 
where Nq is the total number of QSOs used in the anal- 
ysis (Shanks & Boyle 1994; Croom & Shanks 1996). We 
also derive field-to-field errors and compare these to the er- 
rors found in simulations. On small scales, <2 h" 1 Mpc, the 
number of QSO-QSO pairs can be <J10. In this case sim- 
ple root-n errors (Eq. 2) do not give the correct upper and 
lower confidence limits for a Poisson distribution. We use 
the formulae of Gehrels (1986) to estimate the Poisson con- 
fidence intervals for one-sided 84% upper and lower bounds 
(corresponding to lcr for Gaussian statistics). These errors 
are applied to our data for QQ(s) < 20. Above this num- 
ber of pairs root-n errors adequately describe the Poisson 
distribution. 

In our analysis below we will also use the integrated 
correlation function out to some pre-determined radius as a 
measure of clustering amplitude. This is commonly denoted 
by £, where 



£(x)x dx. 



(3) 



As in Paper II we will generally take s max = 20 h^ 1 Mpc 
as this is on a large enough scale that linear theory should 
apply. The effect of z-space distortions due to small-scale 
peculiar velocities or redshift errors is also minimal on this 
scale. 



2.3 Selection functions and incompleteness 

The area of the survey is covered by a mosaic of 2dF point- 
ings. These pointings overlap in order to obtain near com- 
plete coverage in all areas, including regions of high galaxy 
and QSO density. In order to take into account the variable 
completeness between 2dF pointings, due to variations in 
observational conditions, we use a mask that specifies the 
completeness of each survey sector, where we define a sector 
as the unique intersection of a number of circular 2dF fields. 
These masks are fully discussed in Paper XII. The complete- 
ness of each survey strip as a function of angular position 
on the sky is shown in Fig. 2. The distribution of random 
points used in our correlation analysis is constructed to have 
an identical distribution on the sky. In order to minimize the 
influence of low completeness fields, we restrict the analysis 
in this paper to sectors for which the spectroscopic com- 
pleteness is at least 70 per cent. This results in a sample of 
20686 QSOs in the redshift range 0.3 < z < 2.9. 

It is possible that on scales smaller than a 2dF field sys- 
tematic variations in completeness may exist (e.g. see Paper 
XII). In order to test the consequence of these, detailed sim- 
ulations have been carried out (see below) . On larger scales 
small residual calibration errors in the relative magnitude 
zero-points of the UKST plates could add spurious struc- 
ture. These are also assessed using simulations. 

After generating random points according to the an- 
gular distribution specified by the completeness masks, we 
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Figure 1. The distribution of 2QZ QSOs from the final catalogue. The SGP strip is on the left, the equatorial strip on the right. The 
rectangular regions show the distributions projected onto the sky. An EdS cosmology is assumed in calculating the comoving distances 
to each QSO. 
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Figure 2. The completeness map of the 2QZ catalogue for the equatorial (top) and SGP (bottom) regions. The grey-scale indicates the 
percentage of all 2QZ targets that were both observed and positively identified (quality 1) over the two survey strips. 



then assign a random redshift to each point. This random 
redshift is draw from a distribution defined by a polynomial 
fit to the observed n(z) distribution (see Fig. 3a and Section 
3.2.1 below). 

As a direct test of the effectiveness of the above correc- 
tions, we also use random distributions generated by taking 
right ascensions (RAs) and declinations (Decs.) from the 
QSO catalogue itself. We then assign a redshift based on ei- 
ther the fitted n(z) (as above; this we call the RA-Dec mix- 
ing method) or by assigning a random QSO redshift taken 
from the catalogue (the RA-Dec-z mixing method). These 



methods will mimic the 2QZ QSO angular distributions ex- 
actly, but with the effect of reducing the amount of struc- 
ture measured (particularly on larger scales). We examine 
the reduction in large-scale power that these estimates cause 
below. 

These two alternative methods also demonstrate that 
the QSO correlation function is not affected by the deficit of 
close (< 1') pairs in the 2QZ. The deficit is due to the fact 
that the 2dF instrument cannot position two fibres closer 
than ~ 30". It has in large part been alleviated by the over- 
lapping field arrangement in the 2QZ strips, and the fact 
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that the vast majority of QSO pairs which are close in angu- 
lar position have very different redshifts. We therefore make 
no further corrections for this effect in our analysis. 

Extinction by galactic dust will also imprint a signal on 
the angular distribution of the QSOs. Primarily this changes 
the effective magnitude limit in bj by A bj — 4.035xE(B~V) 
where we use the dust reddening E(B — V) as a function of 
position calculated by Schlegel, Finkbeiner & Davis (f998). 
We then weight the random distribution according to the 
reduction in number density caused by the extinction such 
that 

W«t («,«) = 10" -"^j <<*•*>, (4) 

where /3 is the slope of the QSO number counts at the mag- 
nitude limit of the survey. At bj = 20.85, the magnitude 
limit of the 2QZ, the QSO number counts are flat, with 
(3 ~ 0.3. Applying this correction we find that it only makes 
a significant difference to £(s) on scales of ~ fOOO h^ 1 Mpc. 

2.4 Making model comparisons to £(s) 

Below we make comparisons of the data to a number of 
models, both simple functional forms (power laws) and more 
complex, physically motivated, models (e.g. cold dark mat- 
ter; CDM). We use the maximum likelihood method to de- 
termine the best fit parameters. The likelihood estimator is 
based on the Poisson probability distribution function, so 
that 

i=l 

is the likelihood, where v is the observed number of QSO- 
QSO pairs, fi is the expectation value for a given model and 
N is the number of bins fitted. We fit the data with bins 
Alog(s) = O.f, although we note that varying the bin size 
by a factor of two makes no noticeable difference to the resul- 
tant fit. In practice we minimize the function S = — 21n(L), 
and determine the errors from the distribution of AS, where 
AS is assumed to be distributed as x 2 ■ This procedure does 
not give us an absolute measurement of the goodness-of-fit 
for a particular model. We therefore also derive a value of 
X for each model fit in order to confirm that it is a reason- 
able description of the data. In particular this is appropri- 
ate when fitting on moderate to large scales (J>5 h^ 1 Mpc), 
where the pair counts are large enough that the Poisson er- 
rors are well described by Gaussian statistics. 



3 CORRELATION FUNCTION TESTS USING 
MOCK QSO CATALOGUES 

3.1 Mock QSO catalogues 

To test both our correlation function estimation methods 
and the effect of incompleteness we apply our analysis to 
mock QSO catalogues produced from the large Hubble Vol- 
ume simulations of the Virgo Consortium (Frenk et al. 2000; 
Evrard et al. 2002). In particular we make use of the ACDM 
Hubble Volume simulation where data on each particle has 
been output along the observer's past light cone to mimic the 
2QZ. The simulation contains I0 9 particles in a cube that 
is 3000 h^ 1 Mpc on a side. The cosmological parameters 



of the simulation are Ob = 0.04, Ocdm = 0.26, Oa = 0.7, 
H = 70 km s _1 Mpc -1 and cr 8 = 0.9 (at z = 0). The light 
cone data was output in a 75° x 15° wedge oriented along 
the maximal diagonal of the cube, allowing the light cone to 
extend to a scale of ~ 5000 h^ 1 Mpc (z ~ 4). These three 
slices are then split up into 3 largely independent 75° x 5° 
slices, each one mimicing a single 2QZ strip. We note that 
there will be some correlation between the largest structures 
in the different simulation strips, however, it was not prac- 
tical to generate simulations large enough to select many 
completely independent volumes. 

In order to create realistic mock QSO catalogues, the 
mass particles are then biased to give a similar clustering 
amplitude to that observed in the 2QZ (based on the results 
of Croom et al. 2001a). The biasing prescription is based on 
that of Cole et al. (1998) (their model 2), but varying the 
parameters as a function of redshift to match the Croom 
et al. (2001a) results and using a cell size of 20 h' 1 Mpc 
to determine the local density (Hoyle 2000). In our analysis 
below we consider mock catalogues with large numbers of 
biased particles (~ 100000), almost a factor of 10 more than 
a single real 2QZ strip. This allows us to test for possible 
weak systematic affects. Full details of the Hubble Volume 
simulation are given by Hoyle (2000). 



3.2 The effect of different correlation function 
estimates 

There are several issues involved with accurately determin- 
ing the two-point correlation function. We will investigate 
each of these in turn. 



3.2.1 Estimates of the QSO n(z) 

The redshift distributions, n(z), of the two 2QZ slices are 
shown in Fig. 3a. In order to directly compare the two, we 
renormalize the NGP n(z) to contain the same total num- 
ber as the SGP. The two strips have the same overall shape, 
however the we note that they appear to have more struc- 
ture that the n(z) distributions of the Hubble Volume sim- 
ulations shown in Fig. 3b (note that the simulations have a 
cut off imposed at z = 2.2). By examining the spatial dis- 
tribution of the QSOs it is possible to see that the extra 
structure in the n(z) is due to a number of weak large-scale 
structures. For example, the narrow peak in the NGP n(z) 
at z = 1.5 is due to a wall-like feature (top right of Fig. 
1). We must therefore be careful not to remove any excess 
large-scale power by fitting the n(z) on too fine a scale. A 
detailed discussion of structure on very large scales is given 
by Miller et al. (2004). In Fig. 3a we plot the polynomial fit 
(12th order) to the QSO n(z) distribution used to generate 
the random distributions. Tests using higher and lower order 
polynomial fits (8th - 16th order) showed no significant dif- 
ferences between the resultant estimates. We also found 
that different methods of fitting the n(z) of the simulations 
(e.g. spline vs. polynomial) only caused differences at the 
~ 0.1 per cent level, much smaller than the random errors 
in the measurements of £(s) from the 2QZ. 
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Figure 3. QSO and simulation n(z) distributions, a) The n(z) distributions in the two 2QZ slices, SGP (solid line) and NGP (dotted 
line). The NGP has been rcnormalized to the number of QSOs in the SGP to aid comparison. Also shown is the 12th order polynomial 
fit to the combined n(z) (dashed line), b) The n(z) distribution of two Hubble Volume simulation slices each containing 12500 particles. 




s (h~'Mpc) s (h _, Mpc) 

Figure 4. Comparison of masking (filled points) and RA-Dec mixing (open points) methods for the Hubble Volume simulations. Beneath 
each plot we show the ratio of the two correlation function measures, £(s) ma sk/£(s)mixing- a) f (s) measured over a broad redshift range, 
z = 0.3 — 2.2. There is no significant difference between the two estimates, b) £(s) measured over a narrow redshift range, z = 1.35 — 1.70. 
In this case the RA-Dec mixing method produces a correlation function which is ~ 10 — 20 per cent lower than the masking method. 



3.2.2 Masks vs. randomizing 

We next investigate differences between the methods de- 
scribed above to produce the random distributions. In par- 
ticular, although the RA-Dec and RA-Dec-z mixing meth- 
ods are effective at removing any variations in completeness, 
we also need to assess whether they also remove significant 
amounts of large-scale structure. To do this we determine 
the clustering in our simulations using these different meth- 



ods. In Fig. 4 we show a comparison of the masking and 
RA-Dec mixing methods for a single Hubble Volume simu- 
lation slice. When the redshift range is broad (Fig. 4a) there 
is no significant difference between the two methods and the 
ratio of the two (bottom of Fig. 4a) is consistent with 1 at 
all scales. However if we take a narrower redshift interval, as 
in Fig. 4b, we do see significant depression of the clustering 
strength in the RA-Dec mixing method. This is because in a 
narrow redshift interval, the angular clustering of QSOs will 
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Figure 5. Comparison of simulated correlation functions with (open points) and without (filled points) zero-point errors for a) the 
full redshift range and b) a narrow redshift range with z = 1.35 — 1.70. The ratio of the points with and without zero-point errors, 
£(s, <r zp = 0.05)/£(s, cr zp = 0.000), is shown below each plot. 
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Figure 6. Comparison of simulated correlation functions with (open points) and without (filled points) radially dependent incompleteness 
within 2dF fields for a) the full redshift range and b) a narrow redshift range with z = 1.35 — 1.70. The ratio of the points with and 
without radial dependent incompleteness, £(s) r ad/£( s )i is shown below each plot. 



be greater, due to the reduced amount of projection. There- 
fore we conclude that while the RA-Dec mixing method is 
a useful check of the clustering amplitude averaged over the 
full survey, it is not an accurate estimate when measuring 
QSO clustering evolution in narrow redshift slices. The same 
results were found for the RA-Dec-z mixing method. 



3.3 The effect of the survey selection function and 
incompleteness 

We now assess the effect of errors in the survey selection 
function on our estimates of £(s). All these tests are carried 
out using the masking method. Errors in the zero-points of 
the UKST photographic plates are a possible source of excess 
large-scale power. To mimic this effect we divide the simu- 
lated survey strips into 15 5° x 5° regions and apply to each 
a Gaussian random zero-point error Am, with a a — 0.05 
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Figure 7. The QSO £(s) from the 2QZ using the masking method 
(filled points) and and RA-Dcc mixing method (open points). A 
WMAP/2dF cosmology is assumed. Below we show the ratio of 

the two, £(s)mixing/£(s)mask- 



Figure 8. The QSO £(s) from the 2QZ, plotting the results from 
the SGP (filled points) and NGP (open points) separately. A 
WMAP/2dF cosmology is assumed. Below we show the ratio of 
the two, £(s)ngp/£( s )sgp- Note that the scale of the ratio plot 
is broader than the previous similar plots. 



mag. We then modulate the density of sources in that region 
by a factor of i(r a3Am , as the faint end slope of the QSO 
number counts is ~ 0.3. This equates to an error in the QSO 
density of 7 per cent for a zero-point error of 0.1 mag. With 
a — 0.05 the full range of zero-point errors used was ~ 0.15 
mag. We do not expect there to be real zero-point errors 
in the survey larger than this. A comparison of simulated 
correlation functions with and without zero-point errors is 
shown in Fig. 5. We see no systematic differences caused by 
the zero-point errors in either the full redshift interval (Fig. 
5a), or narrower redshift intervals (Fig. 5b). We note that 
if the zero-point errors are increased (to values greater than 
the likely photometric errors in the survey) then significant 
differences can be seen. With a — 0.1 mag there are system- 
atic offsets in £(s) at the level of ~ 1 per cent which become 
significant on scales greater than ~ 40 h~ x Mpc. 

Another possible cause of systematic errors in £ (s) is the 
variations in completeness within 2dF fields. These can be 
caused by systematic errors in astrometry or field rotation 
which will be worse at the edges of a field, or atmospheric 
refraction effects, if a field was observed at a different hour 
angle to that which it was configured for. Paper XII showed 
that although radially dependent completeness is noticeable 
when observations of many individual fields are averaged to- 
gether, if the overlap between fields and repeat observations 
are taken into account there is no systematic decline in com- 
pleteness towards the edge of 2dF fields. In order to confirm 
that completeness variations within 2dF fields will not im- 
pact on our clustering analysis we perform detailed tests. 
We first position our 2dF field centres along the simulation 
strips, and then apply spectroscopic completenesses selected 
randomly from the actual field completenesses found in the 
survey. A mask is also generated to correct for this variable 
incompleteness. We then modulate the completeness within 
each simulated 2dF field such that it mimics the radial de- 



crease seen in Paper XII (filled points in their Fig. 18). We 
then calculated £(s) from these simulations, using a com- 
pleteness mask which corrects for all effects apart from the 
variation in completeness within the 2dF fields. This is a 
worst case scenario, as in the simulations we allocate an ob- 
ject to only one field, and then derive the radial completeness 
variation from the centre of that field. In the actual survey, 
objects without IDs could be observed in overlapping fields. 
We compare the results to £(s) measured without the radial 
completeness variations in Fig. 6. We find that the radial 
completeness variations have no significant impact on £(s) 
for either the whole redshift range or in narrower redshift 
intervals. We also determine the effect of radial incomplete- 
ness on £(s) in narrow redshift intervals (which is used ex- 
tensively in Section 5). The radial incompleteness typically 
only changes £(s) by 2 — 5 per cent, with the worst case 
being 10 per cent. Given that the radial selection model is a 
worst case scenario, and that the measurement errors in £(s) 
are at least 20 per cent, any radial dependence of complete- 
ness within 2dF fields will not impact on our conclusions 
presented below. 



4 THE REDSHIFT AVERAGED QSO 
CORRELATION FUNCTION 

The above simulations confirm that our methods of corre- 
lation analysis, and any residual systematic errors in the 
2QZ should not significantly bias our estimates of £(s). We 
now present the results of applying our correlation analysis 
to the final 2QZ sample, beginning with £(s) averaged over 
the redshift range 0.3 < z < 2.2, for the most part, assum- 
ing a WMAP/2dF cosmology. We note that here we restrict 
the redshift range to regions of high completeness, and do 



© 0000 RAS, MNRAS 000, 000-000 



The 2QZ two-point correlation function 9 



0.03 
0.02 

0.01 



-0.01 



n m = 0.27 n A =0.73 
0.3<z<2.2 



1 



r 



i 



500 



1000 
s (h~ 1 Mpc) 



I 500 



2000 



Figure 9. The QSO £(s) from the 2QZ on scales — 2000 h 1 Mpc, plotted on a linear scale. Error bars are derived from the field-to-field 
variance between six sub-samples of the data set. 



not include QSOs above z = 2.2. This is because the mean 
QSO colours move progressively further into the stellar locus 
above this redshift making the sample increasingly sensitive 
to small systematic errors in selection. This sample contains 
18066 QSOs and has a mean redshift of z — 1.35. 



4.1 Results 

We first plot a comparison between the masking method and 
the RA-Dec mixing method for the redshift averaged QSO 
£(s). This is shown in Fig. 7. Note that we only plot £(s) on 
scales greater than 1 ft" 1 Mpc as we find no QSO-QSO pairs 
on scales smaller than this (in a WMAP/2dF cosmology). 
Also, for any other bins without QSO-QSO pairs we plot 
a point on the bottom x-axis without an error bar. We see 
that on all scales the two estimates are consistent within the 
Poisson measurement errors. There is some indication that 
the RA-Dec mixing method is slightly systematically lower 
than the mask method on scales > 20 h~ x Mpc, which could 
be an indication of a weak systematic error in the mask 
method, but this is not a significant deviation. Given the 
consistency of the two methods, unless we state so explicitly, 
we will use the mask method for all of our £(s) estimates. 

In a second check of the consistency of our results we 
plot a comparison of the measured £(s) in each of the NOP 
and SGP strips (Fig. 8). Although, the £(s) measured from 
the two strips is in broad agreement, the NOP strip shows 
slightly stronger clustering on scales > 20 hT 1 Mpc. Com- 
paring the estimates of £(s) on different scales in the two 
strips we find that they are consistent (0.9u, 1.3u and 0.6cr 
differences for s = 20, 30 and 50 hT 1 Mpc respectively). 

The large volume probed by the 2QZ allows £(s) to be 
probed on very large scales, in excess of ~ 1000 h- 1 Mpc. 
Most models do not predict any signal in £(s) at large scales, 
however, there have been some claims of features in the QSO 
£(s) (including using data from the 2QZ). E.g. Roukema, 
Mamon & Bajtlik (2002) claimed to see several features, in- 
cluding a positive feature at the level of ~ 4 per cent on a 
scale of ~ 240 h' 1 Mpc in the £(s) of ~ 2300 QSOs from 
the initial release of 2QZ catalogue (Croom et al. 2001b). 
To test these claims we make an estimate of the 2QZ £(s) 



to the maximum scales probed by the sample. The results 
of this are shown in Fig. 9 for the WMAP/2dF cosmology 
(Roukema et al. assume Q m = 0.3 and Qa = 0.7, but our re- 
sults are similar for both cosmologies). As Fig. 9 probes very 
large scales, where QSO pairs could be correlated, we deter- 
mine errors by measuring the variance between six subs- 
regions of the full data set (three 5° x 25° regions in each 
2QZ strip) . The errors plotted is the measured rms between 
the six subsamples divided by \/6 to account for the greater 
volume of the full sample. We note that on the largest scales 
even these field-to-field errors will somewhat inaccurate. By 
comparing the QSO-QSO pair counts for the full region and 
the six sub-regions we find that at ~ 200 h^ 1 Mpc ~ 10 per 
cent of pairs come from correlations between different sub- 
regions. By ~ 1000 hT 1 Mpc this number has risen so that 
approximately half of all QSO-QSO pairs are from QSOs in 
different sub-regions. This means that on large scales there 
will be significant correlation between the sub-regions, but 
the reduction of pairs in each sub-regions will also increase 
the Poisson noise. 



There is little evidence of any strong deviation from zero 
on any scale larger that ~ 100 h~ Mpc and the QSO £(s) 
is zero to within 0.5 per cent over a broad range of scales. 
One point (at 90 /i -1 Mpc) deviates from zero by ~ 1 per 
cent. There is no evidence for a feature at ~ 240 h~ l Mpc. At 
various different scales there are some points that are greater 
than la from zero. A x 2 test comparing the data to £(s) = 
at s = 100 - 1000 h' 1 Mpc gives x 2 = 76.1 with 45 degrees 
of freedom (dof), which implies significant deviations at the 
99.7 per cent level. The rms scatter over this scale range 
is ±0.002. The level of deviations away from zero at large 
scales is so small that we cannot be confident that they are 
real features and not due to low level residual systematics. 
However, residual systematic effects at this level will not 
affect any of our conclusions and we can have confidence that 
the masks used to define the selection function are removing 
structure not due to QSO clustering. 
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Table 1. The results of power law fits to the 2QZ averaged over the redshift range 0.3 < z < 2.2. Model fits assuming a power 
law in 2-space [(s/so) -7 ] and a power law in real-space [(r/rrj) '] are presented (the second for a WMAP/2dF cosmology only). The 
real-space power law is corrected for the effects of linear and non-linear 2-space distortion. We list the cosmology assumed, the scales fit 
over, the best fit parameters and associated errors, the measured x 2 values, number of dof, u and probability of acceptance, P(< \ 2 )- 



Model n m ,f2 A Smin.Smax s /r 7 X 2 v P(<X 2 ) 

(s/s )-' 
(s/so)-> 

(r/r )-< 
(r/ro)-< 



0.27,0.73 


1.0,100.0 


5.55 


0.27,0.73 


1.0,25.0 


5.48 


1.00,0.00 


1.0,100.0 


3.89 


1.00,0.00 


1.0,10.0 


3.88 


0.27,0.73 


1.0,100.0 


5.81 


0.27,0.73 


1.0,25.0 


5.84 



+0.29 

-0.29 
+0.42 
-0.48 
+0.18 
-0.18 
+0.43 
-0.53 
+0.29 
-0.29 
3 



1.633 
1.20 



+0.054 
0.054 
+0.10 
0.10 



1 71 o+0.052 
± -' ±J -0.052 

0.86lg;l? 
1.647+°-^ 



37.7 
8.1 
42.6 
5.6 
20.4 



18 
12 
18 
8 
18 
12 



4.6c-3 
7.8c-l 
9.2c-4 
7.0c-l 
3.1e-l 
8.4e-l 



4.2 Fitting models to the QSO £(s) 

We now attempt to fit a variety of models to the data. The 
simplest model traditionally fitted to correlation function 
estimates is a power law of the form 

= (6) 

where so is the comoving correlation length, in units of 
h^ 1 Mpc. We first fit a power law over the full range of 
scales where significant clustering is detected, from 1 to 
100 h- 1 Mpc, using the maximum likelihood technique. For 
the WMAP/2dF cosmology, this resulted in best fit param- 
eters (s , 7) = (5.55 ± 0.29, 1.633 ± 0.054), however this fit 
is unacceptable at the 99.5 per cent level (see Table 1). This 
best fit power law (solid line) is compared to the data in 
Fig. 10a and it can be seen that the data are flatter on small 
scales and steeper on large scales than model. We then vary 
the maximum scale that we fit. Only by reducing this to 
~ 25 h- 1 Mpc is an acceptable power law fit achieved. Over 
the range 1 — 25 h^ 1 Mpc we find best fit values (so, 7) = 



(5A8t ^, 1.20+01,,). The power law slope is significantly 
flatter when the fit is performed on these smaller scales, but 
the scale length, so is largely unaffected. This shows that the 
shape of the QSO £(s) changes with scale and does not fol- 
low a single pure power law, but steepens at large scales. We 
also fit similar power law models to £(s) estimated assum- 
ing an EdS cosmology. Over the range s = 1 — 100 h^ 1 Mpc 
we find (s , 7) = (3.89 ± 0.18, 1.713 ± 0.052), but as for the 
WMAP/2dF cosmology, this is clearly rejected (at the 99.9 
per cent level) (see Fig. 10b). As above, fitting on a more 
restricted range of scales allows acceptable fits. We find an 
acceptable power law fit on scales s = 1 — 10 h^ 1 Mpc with 
(s , 7) = (3.88^^,0.8612:1?) (see Fig. 10b). The apparent 
break in the QSO £(s) is unsurprising given that we gener- 
ally only expect power law clustering in the regime where 
clustering is non-linear. Similar breaks have been seen in the 
clustering of low redshift galaxies (e.g. Hawkins et al. 2003). 
On scales J>10 hT 1 Mpc where £(s) < 1 clustering should 
be close to linear. Other affects, such as z-space distortions 
could also distort the measured £ (s) away from a power law. 
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Figure 11. a) The difference between redshift measurements for 
repeated QSO observations in the 2QZ (using only quality 1 iden- 
tifications and redshifts) as a function of mean redshift (points), 
also shown is the calculated rms in Az = 0.1 bins (solid lines), 
b) The rms redshift difference divided by 1 + z as a function of 
mean redshift. The mean (7 Z /(1 + z) is shown by the dashed line. 



We assess the impact of 2-space distortions on a power 
law. Small scale peculiar velocities will tend to reduce £(s) 
on small scales. Both intrinsic peculiar velocities and redshift 
measurement errors will generate a similar effect. If due to 
intrinsic peculiar velocities, this should be best described by 
an exponential distribution (Ratcliffe et al. 1998; Hoyle et 
al. 2002; Hawkins et al. 2003) such that 



/exp {w z 



1 



1/2 



exp 



-V2 



(7) 



where (w,) 1 ' 2 is the rms pairwise line-of-sight velocity dis- 
persion. If it is the redshift measurement errors which dom- 
inate, then the distribution may be better described by a 
Gaussian, 



/normftfz) — 



exp - 



2{wi) 



(8) 



The rms pairwise redshift error measured from repeat ob- 
servations of 2QZ QSOs is given as a z = 0.0027z in Paper 
XII. We have re-assessed this redshift error using the same 
data as Paper XII (Fig. 11) and find that a better esti- 
mate of the pairwise redshift error is a z = 0.0014(1 + z) 
(the dashed line in Fig. lib). Thus the pairwise velocity 
error [Sv = cSz/(l + z)] corresponding to this redshift er- 
ror is 5v z = 416 km s _1 largely independent of redshift. To 
this we need to add the intrinsic velocity dispersion of the 
QSOs, 5vi. At low redshift the typical intrinsic galaxy pair- 
wise velocity dispersion is ~ 500 km s _1 (e.g. Hawkins et 
al. 2003) at z ~ 0.15. We note that Hawkins et al. did not 
include the factor of 1 + z in Eq. 11 (see below). Correcting 
for this, the pairwise velocity is actually ~ 430 km s _1 . It 
is uncertain whether this will decline with redshift. While 
the dark matter velocity dispersion should decline, as QSOs 
are biased tracers of large-scale structure, their pairwise ve- 
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Figure 12. a) model power law correlations functions demon- 
strating the effects of z-space distortions, assuming a real space 
?( r ) = (f/5)~ 1,8 . In each case £ and £ are plotted with £ being 
the upper line. We show £(r) (solid line), £ii n (s) (dotted line) and 
and £non-lin( s ) (dashed line). For the redshift-space distortion 
model we assume a WMAP/2dF cosmology, at a mean redshift 
of 1.35 with P(z) = 0.4 and (w?) 1 / 2 = 690 km s _1 . b) the ra- 
tio of different models comparing the ratios of §non-lin( s )/£( r ) 
(upper dashed line) and f n0 n— lin(s)/£ii n (s) (upper dotted line). 
The other two dashed and dotted lines are the £ equivalents. The 
two solid lines are set at 1.0 and at (1 + 2/3/3 + /3 2 /5) = 1.30 for 
P = 0.4. 



locity may not decline. Zhao, Jing & Borner (2002) pre- 
dict that the pairwise velocity dispersion of Lyman-break 
galaxies at z ~ 3 could be ~ 200 — 400 km s _1 . Given 
the uncertainty in the evolution of Svi we will assume a 
fixed value of ~ 430 km s _1 at all redshifts, noting that 
any evolution is likely to reduce this value. A final issue 
that needs to be considered is the velocity error due to in- 
trinsic emission-line shifts in QSOs, 5v\. The UV emission 
lines in QSO spectra typically show blue-shifts relative to 
their systemic velocity, this is particular so of lines such as 
CIV. Richards et al. (2002) demonstrated that the dispersion 
between centroids of CIV and Mgll lines was 511 km s _1 , 
while the dispersion between Mgll and [OIII] was a some- 
what smaller 269 km s _1 . This dispersion will cause an ex- 
tra dispersion in our redshift estimates which is not taken 
into account by the repeat observations (as they are repeats 
of the same QSO spectrum). Thus 8v\ should take values 
in the range 200 — 450 km s -1 allowing for measurement 
errors (Richards et al. 2002). Combining the three compo- 
nents of velocity dispersion together in quadrature results in 
(w 2 ) 1//2 ~ 630 — 750 km s _1 . In our analysis below we will 
assume a value of 690 km s _1 which lies in the middle of 
this range. As a combination of 5v\ and Svz dominates the 
total pairwise velocity dispersion, we use Eq. 8 to model the 
effects of z-space distortions on small scales. We note that 
other authors (e.g. Outram et al. 2004; Hoyle et al. 2002) 
used a similar value of (w 2 ) 1 ^ 2 ~ 800 km s _1 (however they 
miss the factor of 1 + z in Eq. 11 below). 
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We should also take into account the effect of linear 
z-space distortions. Kaiser (1987) showed that 

£( S )^(r)(l + ^+i/3 2 ), (9) 

where £(r) is the real-space correlation function and j3 — 
fl^'/b. More generally, £(cr, 7r), the correlation function 
across (the a direction) and along (the tt direction) the line 
of sight is distorted, such that 

l o — 7 

assuming that £(r) is a power law (Matsubara & Suto 1996). 
fx is the cosine of the angle between r and tt (the distance 
along the line of sight), and 7 is slope of the power law. 
Then including the effects of non-linear z-space distortions, 
the full model for £ (a, it) is given by 

/oo 
£'[cr, 7T - (1 + z)w z /H(z)]f noTln (w z )dw z , (11) 
-oc 

where £'[<t,it — (1 + z)w z /H(z)] is given by Eq. 10, /„ orm (m z ) 
is given by Eq. 8 and H(z) is Hubble's constant at a red- 
shift, z. Finally, we carry out a spherical integral over the 
model £(<r, it) to derive the model £(s) which we then fit 
to the data. We note that there is an extra factor of 1 + z 
in Eq. 11 compared to previous works (e.g. Hawkins et al. 
2003; Hoyle et al. 2002). This is because the velocity disper- 
sions are generally given in proper coordinates, rather than 
comoving coordinates. At low redshift this has a minimal 
affect, however, at high redshift this extra term boosts the 
effective scale corresponding to a given proper velocity by 
1 + 2 (in fact it approximately cancels out the increase of 
H(z) with redshift, so that the proper velocity dispersion 
corresponds to a similar comoving scale at every redshift). 
It is therefore critical to incorporate this term. In this pa- 
per we are not specifically focussing on £(cr, tt) and z-space 
distortions, but only wish to determine their affect in shap- 
ing the measured £(s). Detailed investigation of £(cr, tt) is 
discussed by da Angela et al. (in preparation). 

Estimates of the strength of z-space distortions via the 
QSO power spectrum have been made by Outram et al. 
(2004). They find that at z = 1.4, the mean redshift of 
the sample used, /3 = 0.4 ± 0.1. We assume this value for 
P and a small-scale velocity dispersion of 690 km s _1 . We 
then produce a grid of model real-space correlation functions 
which are adjusted for these z-space distortions and fitted to 
our observed £(s) using the maximum likelihood technique. 

In Fig. 12a we show a comparison of models with and 
without z-space distortions. Assuming a real-space corre- 
lation function of £(r) = (r/5) -1 ' 8 in a WMAP/2dF cos- 
mology, and the above values of /3 = 0.4 and (w 2 ) 1 ^ 2 = 
690 kms -1 . The solid lines show the real-space £(r) and 
f(r) (see Eq. 3). The model f(r) is a factor of 3/(3-7) = 2.5 
above £(r). The dotted lines show the model £(s) and £(s) 
for linear z-space distortions only (fn n (s) i.e. (3 = 0.4 and 
(w 2 ) 1 ' 2 — 0.0), while the dashed lines show the full model 
with linear and non-linear z-space distortions (£ n on-iin(s) 
i.e. f3 = 0.4 and (w 2 ) 1/2 = 690 km s _1 ). On scales less than 
10 h' 1 Mpc the non-linear z-space distortions cause a sig- 
nificant suppression of £. In Fig. 12b we plot the ratio of 
these various models. The dashed lines are £ n on-im(s) (top) 
and ^non-iin(s) (bottom) divided by £(r) and £(r) respec- 
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Figure 13. The QSO £(s) from the 2QZ (filled points) com- 
pared to the best fit power law model incorporating the affects of 
linear and non-linear redshift-space distortions. A WMAP/2dF 
cosmology is assumed. The fits are carried out on scales s = 
1 - 100 h- 1 Mpc (solid line) and s = 1 - 25 h' 1 Mpc (dotted 
line). 



tively. The dotted lines are £non-iin(s) (top) and £ n on-iin(s) 
(bottom) divided by £n n (s) and fiin(s) respectively. The 
solid lines are set at 1 and at (1 + 2/3/3 + /3 2 /5) = 1.30 
(for /3 = 0.4). From this it can be seen that on scales 
~ 20 — 30 h~ Mpc and larger the affect of non-linear 
z-space distortion is small, while the linear term affects 
£ on all scales. For the above power law, we find that 
fnon-iin(s)/£ii„(s) = 0.93, 0.97 and 0.99 for s = 20, 30 and 
50 h^ 1 Mpc respectively. 

To begin with we assume a power law model for £(r) 
(Eq. 6). We generate a grid of models with different power 
law slopes (7), and fit these models to the data using the 
maximum likelihood technique over the range s — 1 — 
100 h' 1 Mpc. The resulting best fit model with (3 = 0.4 
and (w 2 ) 1 / 2 = 690 km s _1 is shown by the solid line in Fig. 
13. We find a power law slope of 7 = 1.866 ± 0.060 and 
a real-space scale length ro = 5.81 ± 0.29 h^ 1 Mpc. This 
provides an acceptable fit to the data with \ 2 = 20.4 (18 
dof) and an acceptance probability of 31 per cent. If we fit 
over a more restricted range of scales, noting that we ex- 
pect deviations from a pure power law in real space on large 
scales, then we find best fit values of 7 = 1.647 ± 0.047 and 
r = 5.84 ± 0.33 h' 1 Mpc for s = 1 - 25 h' 1 Mpc. Both 
fits are compared to the data in Fig. 13 (see also Table 1). 
When fitting on smaller scales the power law slope is flatter, 
however, ro is unchanged. It can be seen that the affect of 
small scale z-space distortions has a significant impact on 
scales less than ~ 10 ft -1 Mpc. 

More generally we should fit a model where the shape 
of £(r) is governed by the underlying physics of the dark- 
matter distribution (e.g. CDM). In particular, Hamilton et 
al. (1991,1995) provide an analytic description of the generic 
linear CDM £(r). The input parameters for the CDM model 
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Figure 14. The QSO from the 2QZ (filled points) compared 
to a WMAP/2dF CDM model. The models shown are the linear 
real-space mass correlation function, £p(r) (short dashed line), 
the non-linear £p(r) (dotted line). The non-linear £p(r) scaled 
by the best fit bias value (long dashed line) and the non-linear 
mass correlation function corrected for z-space distortions, £p(s), 
scaled by the best fit bias (solid line). 



are taken from the now standard WMAP/2dF cosmologi- 
cal model (Spergel et al. 2003; Percival et al. 2002) with 
fi m = 0.27, fi A = 0.73, n b = 0.04, H = 71 km s" 1 MpcT 1 , 
as = 0.84 (at z = 0). We calculate the model £(s) at the 
mean redshift of the 2QZ sample (z = 1.35), and correct 
for the affects of non-linear clustering (Hamilton et al. 1991; 
Jain et al. 1995). Linear and non- linear z-space effects are 
accounted for as above, but using the more general prescrip- 
tion of Hamilton (1992) rather than Eq. 10 for the linear dis- 
tortions. For the z-space distortions we assume /3 = 0.4 and 
(w^) 1 ' 2 = 690 km s -1 . We then perform a maximum like- 
lihood fit for a single parameter, a scale independent QSO 
bias, over the scale range s = 1 — 100 h^ 1 Mpc. QSO bias 
is defined as 



Mr) 



(12) 



where £cj(r) and £p( r ) are the real-space QSO and mass 
correlation functions respectively. We note that our assumed 
value of [3 includes an implicit assumption of QSO bias. If 
we substitute the £,q(t) m Eq. 12 with that from Eq. 9 and 
solve the resultant quadratic in 6q(«) we find that 



fo(«) 4^ 2 (z) 

Mr) 45 3 



(13) 



This relation thus directly gives us the QSO bias at a red- 
shift z, but is only strictly true if non-linear z-space distor- 
tions, which affect the shape of £(s), are not present. The 
linear distortions do not affect the shape of £(s) (this is 
exactly the case when there are no non-linear effects, and 
correct to first order in the presence of non- linear effects), 
so we fit a model £ p (s) divided by (1 + 2(3/3 + /3 2 /5) (us- 



Figure 15. The QSO £(s) from the 2QZ (filled points) compared 
to the the 2dFGRS £(s) of Hawkins ct al. (2003) (solid line, with 
±lrr errors shown by the dotted lines). 



ing the same (5 = 0.4 value used above) to obtain the ratio 
%Q(s)/£p(r) seen in Eq. 13. Assuming f2 m (z = 0) = 0.27 
[implying O m (z = 1.35) = 0.83] we find a best fit QSO bias 
of b Q (z = 1.35) = 2.02 ± 0.07. This model is fully consistent 
with the data, with a x 2 — 14.3 from 19 dof (acceptable 
at the 76 per cent level, see the solid line in Fig. 14). The 
implied values of j3 for this best fit bias is j3 = 0.44 ± 0.02. 
This is close to our assumed value of /? = 0.4 and within the 
errors estimated by Outram et al. (2004) of ±0.1. To test 
the impact of making the z-space corrections to our model, 
we also fit the non-linear real space model to the data. This 
results in a best fit bias of 2.12 ± 0.09 (long dashed line in 
Fig. 14), however, this is a slightly worse fit with a\ 2 ~ 25.2 
(19 dof) acceptable at the 15 per cent level. From Fig. 14 we 
see that the real-space model does not have a strong enough 
break at ~ 10 — 20 h^ 1 Mpc to match the data. We con- 
clude that the 2QZ QSO £(s) averaged over redshift is fully 
consistent with the WMAP/2dF cosmology once allowance 
is made for the affects of z-space distortions. 



4.3 Comparisons to other results 

The redshift averaged QSO £(s) from the 2QZ is consistent 
with the current best fit cosmological model, after allowing 
for a linear bias 6q(z = 1.35) = 2.02 ± 0.07. We now com- 
pare our results to those from other estimates of £(s). We 
find that there is very good agreement between the 2QZ £(s) 
and 2dF Galaxy Redshift Survey (2dFGRS; Hawkins et al. 
2003) £(s) both in the shape and amplitude (see Fig. 15). 
We note that the 2QZ £(s) may be slightly flatter than that 
of the 2dFGRS on small scales, as would be expected given 
the smaller influence of non-linear clustering at high redshift 
together with the larger impact of non-linear z-space distor- 
tions. However this is not significant. While the agreement 
in shape is not particularly surprising, the impressive match 
in amplitude is more surprising. This was also found in the 
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Figure 16. The QSO £(s) from the 2QZ (filled points) compared 
to that for NVSS detected 2QZ QSOs (open points). The radio 
detected £(s) uses broader bins of A log(s) = 0.2. 



preliminary 2QZ data release (Croom et al. 2001a). Consid- 
ering the evolution of clustering seen (see Section 5 below) , 
this must be considered as something of a coincidence. 

A number of authors have measured the spatial cluster- 
ing of radio galaxies over a range of redshifts. Overzier et al. 
(2003) finds a real-space clustering scale-length ro = 14 ± 
3 h^ 1 Mpc at z ~ 1 for powerful radio galaxies, while weaker 
radio sources appear less clustered, with ro ~ 4 — 6 h^ 1 Mpc. 
The clustering of 2QZ QSOs (which are largely radio quiet) 
is more similar to the radio weak sources. The 2QZ contains 
a small fraction of sources detected in the radio. There are 
428 2QZ QSOs in the redshift range 0.3 < z < 2.2 that are 
detected by the NRAO VLA Radio Survey (NVSS; Condon 
et al. 1998). The £(s) we measure for this radio-detected pop- 
ulation is shown in Fig. 16 (open circles). The small number 
of sources and their low surface density means that there 
is barely a detection of clustering, with only 2 QSO pairs 
detected vs. 1.15 expected at s < 20 ft -1 Mpc. The cluster- 
ing of radio-detected QSOs in the 2QZ does not therefore 
impact on the clustering measurements of the full sample. 
There is a clear difference between the clustering of radio- 
quiet QSOs, as sampled by the 2QZ, and powerful radio 
galaxies, implying that radio galaxies must exist in more 
massive dark matter halos that radio-quiet QSOs. 

The low redshift galaxy cluster correlation function has 
a much higher amplitude with so typically 12 — 25 h^ 1 Mpc 
depending on the richness of the clusters (Bahcall et al. 
2003). There are few measurements of the cluster correla- 
tion length at high redshift. Gonzalex, Zaritsky & Wechsler 
(2002) find that approximately velocity dispersion limited 
samples of clusters at z = 0.35 — 0.575 have similar cluster- 
ing scale lengths to local clusters. For a WMAP/2dF cos- 
mology, linear theory predicts that the amplitude of mass 
clustering between z = 1.35 and z = will increase by a 
factor of ~ 3.4, which is equivalent to an increase in so by 
a factor of 2.0 (assuming 7 = —1.8). Hence, even if QSO 



clustering at a mean redshift of z = 1.35 evolved as strongly 
as linear theory evolution allows (making no allowance for 
evolution of bias), the descendents of objects that contained 
QSOs at z ~ 1.4 could not be clustered any more strongly 
than poor clusters at low redshift. Below we make a more de- 
tailed analysis of the evolution of QSO clustering to extend 
this analysis. 



5 THE EVOLUTION OF QSO CLUSTERING 

Above we have calculated £(s) average over a broad red- 
shift range. Under the assumption that QSO bias is largely 
scale independent (at least compared to the uncertainties in 
the clustering measurements) this should preserve the cor- 
rect underlying shape of £(s), particularly on large scales. 
However, according to the standard picture of gravitational 
growth of structure, the mass distribution should evolve with 
redshift. Croom et al. (2001a) showed that QSO cluster- 
ing was constant or slightly increasing with redshift, with 
so — 5 h^ 1 Mpc up to z ~ 2.5. This demonstrated that 
QSOs must be biased tracers of the matter distribution, and 
that the amount of bias must evolve with redshift. Below we 
repeat this analysis with the final 2QZ data set, and discuss 
in detail the implications for QSO formation models. We will 
assume a WMAP/2dF cosmology unless stated otherwise. 



5.1 Measurements of £(s, z) 

We split the QSOs up into 10 redshift intervals, such that 
there are approximately equal numbers of QSOs (~ 2000) 
in each bin. Here we sample the redshift range 0.3 < z < 2.9 
and note that the final redshift interval 2 = 2.25 — 2.90 
could be affected by systematic variations in completeness 
on large scales. We perform the correlation analysis as de- 
scribed above on each of these sub-samples. In particular we 
use the mask method to correct for incompleteness, as the 
RA-Dec mixing method was shown to significantly suppress 
clustering measurements in narrow redshift intervals (see 
Section 3.2.2). We do, however, perform tests with the RA- 
Dec and RA-Dec-z mixing methods to confirm that there are 
no obvious unaccounted for systematic errors in our analy- 
sis. The resulting correlation functions are plotted in Fig. 
17. 

In order to make quantitative measure of the clustering 
properties we calculate £ (20) (Eq. 3) for each redshift inter- 
val. To test for any evidence of a change in shape of £(s) 
we also calculate £ using radii of 30 and 50 h' 1 Mpc. The 
evolution of £ is plotted in Fig. 18a using all three scales 
(the values are also listed in Table 2). In each case there is 
a general trend for £ to increase with redshift. To assess the 
significance of the evolution we perform a Spearman rank 
correlation test on the £ values. We find Spearman rank- 
order correlation coefficients, p — 0.721, 0.648 and 0.552 for 
£ determined at a radius of 20, 30 and 50 h~ l Mpc respec- 
tively. These correspond to correlation significances of 98.1, 
95.7 and 90.2 per cent. We note, of course, that as these are 
integral measures they are not independent of each other. 
The above test implies a significant correlation with red- 
shift, however the data are still found to be consistent (via 
a X 



2 test) with a single parameter model which is constant 
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Figure 17. The QSO £(s) from the 2QZ (filled points) as a function of rcdshift in 10 rcdshift bins containing approximately 2000 QSOs 
each. The best fit power law is shown in each case (solid line), as well as the best fit fixing 7 to be 1.20 (dashed lines). We also show the 
best fit power law for the full rcdshift range (0.3 < z < 2.2) for comparison (dotted line). A WMAP/2dF cosmology is assumed. 
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Table 2. 2QZ clustering results as a function of rcdshift for a WMAP/2dF cosmology. All fits are on scales s = 1 — 25 Mpc. We list 
the redshift interval, and mean redshift, apparent magnitude and absolute magnitude (assuming h = 0.71) for each bin together with 
the number of QSOs used. The best fit values of so ( m comoving units of h~ 1 Mpc) and 7 are given with their x 2 values, number of 
dof, v and probability of acceptance, P(< X )■ Lastly we also list the measured values of £(s) for s = 20, 30 and 50 h~ 1 Mpc. 
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Figure 18. a) The evolution of £(s) for three different values of s = 20, 30 and 50 h^ 1 Mpc (open circles, filled circles and open squares 
respectively). There is evidence for an increase in £ with increasing rcdshift in all cases, b) The ratios of £ as a function of rcdshift for 
£(20)/§(30) (open circles) and ^(50)/^(30) (filled circles). The redshift averaged mean values for the ratios arc indicated by the solid 
lines. The ratios are consistent with an unchanging shape for £(s). The WMAP/2dF cosmology is assumed. Also plotted are the expected 
ratios for a CDM model with WMAP/2dF parameters (dotted lines). 



with redshift (only rejected at the 81, 77 and 75 per cent 
levels for £(20), £(30) and £(50) respectively). 

In Fig. 18b we show the ratio of f(20)/f(30) and 
£(50)/£(30) to provide a simple test for any evidence of 
a change in the shape of £(s) with redshift. These ratios 
are consistent with being constant over the full redshift 
range of the data set, suggesting that the shape of £(s) does 
not change significantly with redshift. We also compare the 
£ ratios to those assuming a CDM power spectrum in a 
WMAP/2dF cosmology (dotted lines in Fig. 18b). These 
are fully consistent with the observed ratios. In Fig. 19 we 
show the evolution of £(s) for an EdS cosmology. In this 
cosmology clustering is completely constant as a function of 



redshift, a Spearman rank correlation test shows no signifi- 
cant correlation. 

We next fit a simple power law model (Eq. 6) . In Section 
4.2 we find that a power law is an acceptable fit to the 
redshift averaged QSO £(s) on scales s = 1 — 25 h^ 1 Mpc. 
We therefore fit the data sub-divided into redshift intervals 
over the same range of scales. The best fit so and 7 values 
are shown in Fig. 20 (and listed in Table 2). We carry out 
a Spearman rank test on both so and 7 vs. redshift. For 
so we find p — 0.770 (99 per cent significant), while for 7 
we find p = —0.030 (7 per cent significant). The measured 
values of so are inconsistent with a constant value at 98 per 
cent significance. Given the lack of evolution in 7 we now 



© 0000 RAS, MNRAS 000, 000-000 



The 2QZ two-point correlation function 17 



0.6 




Figure 19. The evolution of for three different values of 
s = 20, 30 and 50 h~ 1 Mpc (open circles, filled circles and open 
squares respectively) in an EdS cosmology. 
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Figure 20. a) The fitted values of so with freely varying 7 (filled 
points) and fixed 7 (open points) as a function of rcdshift. b) The 
best fit values of 7 as a function of rcdshift. The dotted lines indi- 
cate the best fit values to the full redshift range. A WMAP/2dF 
cosmology is assumed. 



fix its value and re-perform the fitting. For this we use the 
best fit power law slope of 7 = 1.20. The so values derived 
are plotted in Fig. 20 (open points). These are similar to 
those found when allowing 7 to vary freely. A Spearman 
rank correlation test confirms that the correlation is still 
present with p — 0.842 significant at the 99.8 per cent level. 

Examining the highest redshift bin in Fig. 17 we see 
that there is significant signal at scales ~ 70— 100 h^ 1 Mpc. 



This redshift interval at 2.25 < 2 < 2.90 has a large vari- 
ation in completeness with redshift, as the absorption due 
to the Lyman-cv forest quickly moves the mean QSO colours 
into the stellar locus (see Paper XII). We do not need to cal- 
culate the absolute completeness in each redshift interval, as 
we rely on fitting to the observed shape of the QSO n(z) re- 
lation. However, if this fit is not accurate enough over a given 
redshift interval, or there are systematic differences in the 
n(z) covering different regions of the 2QZ survey, extra spu- 
rious large-scale structure could be added. We test for the 
presence of any such systematic affect by first calculating the 
£(s,z) using RA-Dec-z mixing. This produces estimates of 
£(s) which are systematically biased low (see Section 3.2.2), 
however any broad trends should still be present. We find 
that the highest redshift bin still has the largest best fit value 
of s using these mixing methods. As a second test we cal- 
culate £(s) for the 2.25 < z < 2.90 interval by normalizing 
the total number and the redshift distribution of the ran- 
dom points within each UKST field. This would remove the 
affects of any UKST photometric zero-points errors or the 
differential affects of variability on completeness in differ- 
ent fields. The results of this analysis are indistinguishable 
from those using masking and the full 2QZ strips. While it is 
possible that this excess large-scale structure is still caused 
by systematic error, its size does not influence any of our 
main results below. Infact the final redshift bin could be 
completely ignored without changing our basic conclusions. 

5.2 Comparison to simple models 

Following Paper II we test a number of simple models 
against the observed data. To be conservative we use the 
£ (20) measurements, rather than the best fit so values which 
are dependent on the range of scales fit and assumptions 
concerning the slope, 7. We note that removing the high- 
est redshift point does not remove the detected correlation 
between £(20) and redshift, although it does reduce its sig- 
nificance (p = 0.617, significant at the 92 per cent level). 
The significance of the correlations of £ (30) and £(50) with 
redshift are also reduced removing when the highest redshift 
point is removed (to 85 and 69 per cent respectively). 

We compare our results to the expected growth in den- 
sity perturbations from linear theory, which should be ap- 
plicable on the scales we are probing. For an EdS universe, 
the linear growth rate, D(z), is given by D(z) — 1/(1 + 2), 
and for other cosmologies we use the accurate fitting formula 
of Carroll, Press & Turner (1992). In Fig. 21a we plot the 
measured £(20) for QSOs vs. linear theory models (dotted 
lines). We assume a CDM model with WMAP/2dF param- 
eters. In this model the values of £(r, 2 = 0) for the mass 
distribution are 0.254, 0.123 and 0.042 for r = 20, 30 and 
50 h^ 1 Mpc respectively. We plot two linear theory lines, 
the first (lower dotted line) assumes the above normaliza- 
tion given by WMAP/2dF, which is significantly below the 
points at all redshifts. The second (upper dotted line) is 
the linear theory model re-normalized by a constant bias to 
a 'best fit' value for the data points. As in Croom et al. 
(2001a) we find linear theory evolution with a fixed bias to 
be in clear disagreement with the data (the probability of 
acceptance is formally 3.6 x 10 -9 ). Assuming an EdS cosmol- 
ogy, we also get a rejection of QSOs following linear theory 
evolution (rejected at the 99.98 per cent level). We next fit 
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Figure 21. a) Our measurement of §(20) for 2QZ QSOs as a func- 
tion of redshift (filled points). The data are compared to linear 
theory gravitational evolution (dotted lines) for two normaliza- 
tions, one normalized to a WMAP/2dF cosmology (lower dotted 
line) and a second normalized to provide a 'best fit' to the data 
points (upper dotted line). We also compare to the best fit for a 
constant §(20) (solid line) and a long-lived model (short dashed 
line), b) The QSO bias, 6q(z) as a function of redshift derived 
from a comparison of §(20) for QSOs to that expected for the 
WMAP/2dF cosmology. The open points are the raw bias val- 
ues [i.e. §Q(s)/§p(r)] while the filled points with error bars are 
the values after making a consistent correction for z-space dis- 
tortions. A simple empirical model is also shown (dotted line), 
c) The mean mass of DMHs containing QSOs derived from the 
measured bias (filled points). We also show the mean mass aver- 
aged over redshift (solid line) and the mean plus twice the rms 
of the points (long dashed line). M*(z), the characteristic mass 
which is just collapsing at a given redshift is denoted by a dotted 
line. The short dashed lines show the median expected growth in 
DMH mass from the mean DMH mass of QSO hosts at z = 0.53, 
1.41 and 2.48. 



the long-lived QSO model discussed by Croom et al. (2001a) 
which has the form 



6q(z) = l + [Mz = 0)-l]/D(z). 



(14) 



This model is equivalent to assuming that QSOs have ages 
of order the Hubble time, and after formation at some ar- 
bitrarily high redshift subsequent evolution is governed by 
their motion within the gravitational potential (Fry 1996). 
It is also equivalent to QSOs forming in density peaks above 
a constant threshold (Croom & Shanks 1996). The best fit 
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Figure 22. The ratio of § non — lin( s )/§lin( s ) (§ including and not 
including non-linear z-space distortions) as a function of redshift 
for (w 2 ) 1 / 2 = 630, 690 and 750 km s" 1 (solid, dotted and dashed 
respectively) at scales s = 20, 30 and 50 h^ 1 Mpc (bottom to 
top). 



value of 6q(z = 0) = 1.64 ± 0.05 (short dashed lines in Fig. 
21a), however, while Croom et al. (2001a) found this model 
was marginally acceptable in a cosmology with fi m = 0.3 
and Qa = 0.7 we find that the extra signal in the final 2QZ 
data set rejects the long-lived model at a significance level 
of 99.97 per cent in the WMAP/2dF cosmology. Fitting this 
model in the EdS Universe gives 6q(z = 0) = 1.40 ± 0.04, 
and is marginally acceptable (rejected at the 89 per cent 
level) . 



5.3 Bias, dark matter halo mass and the evolution 
of QSOs 

By assuming an underlying cosmological model we are able 
to convert the measured values of f to an effective bias by 
making comparisons to linear theory evolution. This allows 
us to directly determine QSO bias as a function of redshift. 
In doing so, we need to account for the affect of z-space 
distortions on the measured values of £(s). The non-linear 
z-space distortions have a small affect on the scales we are 
examining here (see Section 4.2). To determine their affect 
on we derive the ratio of £(s) with linear and non- 
linear z-space distortions to that including only the linear 
distortions, £ n on-iin(s)/£iin(s)- This is plotted for the CDM 
model with WMAP/2dF parameters as a function of redshift 
for {w 2 } 1/2 = 630, 690 and 750 km s" 1 in Fig. 22 (solid, 
dotted and dashed lines respectively). In constructing the 
models we assume values of /3 that are consistent with the 
P(z = 1.4) = 0.4 ± 0.1 of Outram et al. (2004) and also 
account for the evolution of bias we find below. This as- 
sumption of (3 only influences the shape of £(cr, tt) that is 
convolved with Eq. 8 to determine the non-linear z-space 
distortions. Varying the assumed (3 within reasonable limits 
results in negligible difference in the £ n on-iin(s)/£iin(s) ra- 
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Table 3. The derived QSO bias, 6q and DMH mass, M DH as for 
2QZ QSOs at a function of redshift in a WMAP/2dF cosmology. 
We also list the mean redshift and absolute magnitude of each 
redshift interval, as well as the value of M£ derived from the 
polynomial evolution model of Paper XII anc! the space density 
of QSOs, <I>, found by integrating the QSO luminosity function 
between the apparent magnitude limits of the 2QZ. 
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tio (less than 0.5 per cent). We plot the ratio for s = 20, 
30 and 50 hT l Mpc (top to bottom) and see that even at 
s — 20 h^ 1 Mpc the worst correction is only 12 per cent. 
Assuming (w 2 ) 1 ^ 2 = 690 kms -1 , the range of reasonable 
values for (u> 2 } 1,/2 results in a scatter of only ~ 2 per cent 
at s = 20 h^ 1 Mpc and less at larger scales. This is con- 
siderably smaller than the measurement errors in £, and we 
therefore use the derived ratio for (w 2 ) 1 '" 2 = 690 km s _1 
to correct our results for non-linear z-space affects (dotted 
lines in Fig. 22). Linear z-space distortions (Eq. 9) have a 
more significant affect, (e.g. a factor of ~ 1.3 at z ~ 1.4). 
We use Eq. 13 to self-consistently determine the QSO bias 
at a given redshift. 

Fig. 21b shows the derived bias of 2QZ QSOs as a func- 
tion of redshift (filled points) . The open points are the values 
found without accounting for z-space distortions. Here we 
see that QSO bias is strongly evolving with redshift, from 
b Q (z = 0.53) = 1.13±0.18 to b Q (z = 2.48) = 4.24±0.53 (see 
Table 3) . A simple empirical description of the bias evolution 
found is 



b Q (z) = (0.53 ± 0.19) + (0.289 ± 0.035)(1 + zf 



(15) 



which is shown in Fig. 21b (dotted line). At z ~ 0.5 the 
value of 6q is already close to 1, and a simple extrapolation 
of the trend observed would predict that the bias would at 
or below 1 at z = 0. We note at this point that because of 
the apparent magnitude limit of the 2QZ, the mean abso- 
lute magnitude in each interval increases with redshift (see 
Table 2). However, the 2QZ selects QSOs that are close to 
~ Lq (the characteristic luminosity of the QSO optical lu- 
minosity function) at every redshift, and the space density 
of objects in each of the redshift slices is also approximately 
equal. Table 3 lists the values of M bj (z) assuming the poly- 
nomial evolution model of Paper XII (which is an uncertain 
extrapolation beyond z = 2.1). Although the actual values 
of M^(z) should be used with caution as the fitted value 
of M bj (0) is correlated with the bright and faint end slopes 



of the LF, it can be seen that there is little change in the 
relative difference between M^ } (z) and Mb,(z) (less than 1 
mag at z < 2.2). Also listed is the space density found by 
integrating the observed luminosity function over the appar- 
ent magnitude range of the 2QZ for each redshift. Between 
z — 0.5 and z = 2.1 there is only a factor of 2 change in 
space density (increasing to a factor of 2.7 if we include the 
highest redshift bin). Paper XII found that the extrapolated 
M£ (the absolute magnitude equivalent of Lq) at z = has 
in the range —20.5 to —21.6 (where the large range is due to 
correlation between the value of M bj and the bright /faint 
slopes of the QSO LF, and uncertainty in the exact model to 
extrapolate to zero redshift). Thus we would expect that at 
these moderate luminosities, QSOs (or more properly AGN) 
would be close to unbiased at z = 0. It has been shown 
(Hawkins et al. 2003; Verdi et al. 2002) that ~ L* al galaxies 
at low redshift are largely unbiased. This implies that typical 
low redshift AGN (which are much less luminous than those 
at high redshift) are clustered similarly to ~ L* al galaxies. 
There is some direct evidence that this is the Croom 
ct al. (2004c) have shown that the cross-correlation between 
low redshift 2QZ QSOs and 2dFGRS galaxies is equal to the 
auto-correlation of the galaxies. 

Once the bias is derived it is possible to relate this to 
the mean mass of the DMHs that the QSOs reside in. Halos 
of a given mass, M, are expected to be clustered differently 
to the underlying mass distribution. Mo & White (1996) 
developed the formalism for relating mass to bias. This was 
extended to low mass halos by Jing (1998). Both of these 
works were based on the spherical collapse model. Sheth, 
Mo & Tormen (2001) extend the formalism to account for 
ellipsoidal collapse, to provide an improved relation between 
bias and mass. It is this relation that we will use in our 
analysis. The bias is related to the mass via 



b(M,z) 



Va5 c (z) _ 



av 2 yfa + 0.5y / o(a^ 2 ) 1 
{av 2 ) c 



(av 2 Y + 0.5(1 -c)(l-c/2) 



(16) 



where v = S c (z)/a(M, z), a = 0.707 and c = 0.6. S c is the 
critical overdensity for collapse of a homogeneous spherical 
perturbation. For an EdS universe S c = 0.15(12tt) 2/3 ~ 1.69. 
For a general cosmology 8 C has a weak dependence on red- 
shift, which is given by Navarro, Frenk & White (1997). 
a(M) is the rms fluctuation in the linear density field on a 
mass scale, M, and is given by 



k 2 P(k)w 2 (kr)dk, 



(17) 



where P(k) is the power spectrum of density perturbations 
and 



w(kr) 



3(fcrsin(fcr) — cos(fer)) 
(krf ' 



(18) 



which is the Fourier transform of a spherical top-hat of size 

1/3 



- ( 3M 

\ 47rp 



(19) 



po is the mean density of the universe at z = and cor- 
responds to 2.78 x lO n f2 m /i 2 M Mpc" 3 . a(M) at z = 
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is related to that at arbitrary redshift by the linear growth 
factor, D(z), such that 

a(M,z) = a(M)D(z). (20) 

The characteristic mass at any given redshift, M*(z), that 
is, the mass scale which is just collapsing at a given redshift 
is defined by 

'[WW] = jfcy (2D 

We apply Eq. 16 to estimate the typical mass of the 
DMHs containing our QSOs at each redshift. This typical 
mass is plotted in Fig. 21c. We find that the typical Mdh of 
2QZ QSO hosts is largely constant as a function of redshift, 
even though their typical luminosity is increasing at high z. 
There appears a slight tendency for low redshift QSOs to be 
in lower mass DMHs, but a Spearman rank test shows no sig- 
nificant correlation between redshift and Mdh (p = 0.467, 
significant at only the 83 per cent level). The mean mass 
corresponds to Mdh = (3.0 ± 1.6) x 1O 12 /i _1 M (rms er- 
ror). By comparison, the characteristic mass of the Press- 
Schcchter mass function (Press & Schechter 1974), M*, is 
declining quickly at high redshift (dotted line in Fig. 21c). 
M* halos are unbiased (b = 1) at every redshift, with ha- 
los more massive than M* becoming progressively more bi- 
ased. We therefore see that the increasing bias of DMHs 
hosting 2QZ QSOs towards higher redshift makes them in- 
creasingly more massive than M*. However, the increase in 
mass relative to M* is almost exactly cancelled out by the 
evolution of M* to give an approximately constant Mdh. 
We find that Mdh for QSO hosts is in fact very similar 
to M*(z = 0) ~ 3.5 x 10 12 /i _1 M Q . This is effectively the 
same result discussed above, that by extrapolation Lq QSOs 
would be largely unbiased at z ~ 0. The actual mass derived 
is dependent on the exact cosmology used. Varying our as- 
sumed os{z = 0) = 0.84 by ±0.08 [the 2a range from from 
analysis of WMAP and other data (Spergel et al. 2003)] 
gives a range in M DH between (1.52 ± 0.86) x 1O 12 /i _1 M 
and (5.4 ± 2.8) x 1O 12 /i _1 M for a»(z = 0) = 0.76 and 
CTg(z = 0) = 0.92 respectively. Such changes in normaliza- 
tion will affect all redshift intervals equally, and also scale 
the value of M* by an equal amount. So although the de- 
rived mass might be different our overall conclusions (in 
terms of constant Mdh and 6q ~ 1 at z = 0) are not 
affected. Using a different form for the relation between 
6q and Mdh also slightly affects out results. The relations 
described by Mo & White (1996) and Jing (1998) give a 
mean Mdh ^ (1.9 ± 0.9) x 10 12 /i _1 M Q . These show even 
less dependence of Mdh with redshift, as the masses of 
the highest redshift halos are reduced the most. We con- 
firm that similar results are found using the estimates of 
£(30), these give a similar non-evolving Mdh, with a mean 
of (2.2 ± 1.3) x 10 12 h^ 1 Mq. Our mass estimates are consis- 
tent with those derived by Grazian et al. (2004) based on 
the QSO clustering results of Croom et al. (2001a). 

5.3.1 The lifetime of QSOs 

The observation that 2QZ QSOs sample the same mass 
DMHs at every redshift further demonstrates that we cannot 
be seeing a cosmologically long lived population. As the mass 
of DMHs grow with time through the process of accretion 




z 



Figure 23. The 2cr upper limits to QSO lifetime as a function of 
redshift (connected filled circles), based on the growth in mass of 
DMHs. 



and merging, the low redshift descendents of high redshift 
QSOs will inhabit higher mass DMHs, and hence the QSOs 
we observe at high and low redshift cannot be drawn from 
the same single coeval population. We use the formalism for 
DMH evolution developed by Lacey & Cole (1993) to pre- 
dict the median mass of the descendents of DMHs hosting 
QSOs at later epochs. Eq. 2.22 of Lacey & Cole gives the cu- 
mulative probability that a DMH of mass Mi at time t\ will 
merge to form a new DMH of mass greater than M2 by time 
ti. By finding the mass, M2, that corresponds to a probabil- 
ity of 0.5 at a given time ti we have the median mass of de- 
scendent DMHs. In Fig. 21c we plot the evolution of the me- 
dian DMH mass for a starting mass of 3.0 x 1O 12 /i _1 M (the 
mean QSO host Mdh) at z = 0.53, 1.41 and 2.48 (dashed 
lines). At low redshift, there is only limited time for growth, 
and the DMHs of QSO hosts at z ~ 0.5 would only have 
evolved to a mass of ~ 1 x lO 13 /i. _1 M at z = 0. However, 
the highest redshift DMHs hosting QSO have more time to 
evolve and would have typical masses of ~ 6 x 1O 14 /i _1 M 
at z = 0. It therefore appears that 2QZ QSOs at high red- 
shift (z ~ 2) inhabit the progenitors of low redshift galaxy 
clusters, while 2QZ QSOs at lower redshift are located in the 
progenitors of galaxy groups. The growth of Mdh allows us 
to place constraints on the allowable lifetime of QSO activ- 
ity. Low redshift QSOs cannot be the same population of ob- 
jects as at higher redshift if they have masses which are less 
than the mass of the high redshift sources, after accounting 
for their expected growth over time. Therefore calculating 
the time taken to reach the mean QSO host DMH mass plus 
twice the measured rms gives a ~ 2a limit on the lifetime 
of QSO activity (the rms is 1.6 x lO 12 fe _1 M and the long 
dashed line in Fig. 21 shows the mean plus twice this rms). 
The result of this is plotted in Fig. 23 (connected filled cir- 
cles). At high redshift, halos merge more quickly than at low 
redshift, therefore we find that the limits on QSO lifetime 
using this method are smaller at high redshift than at low 
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redshift. At z = 2.48 the 2a upper limit on QSO lifetime is 
6 x 10 8 years, while at redshifts below z = 1.7, the upper 
limit is > 1 x 10 9 years. At z = 0.53 the limit is 3 x 10 9 
years. 

To further constrain QSO lifetimes, a number of authors 
have produced models for QSO clustering in order to try and 
constrain the typical lifetime of QSOs. Martini & Weinberg 
(2001) give fitting functions for their models which relate 
ri, the scale at which the rms fluctuations in the QSO dis- 
tribution is 1 (i.e. aci(r\,z) = 1) to typical QSO lifetime. 
Their model makes some assumptions, including that the 
brightest QSOs are always in the most massive halos at any 
given redshift and that the presence of a black hole is the 
only requirement for QSO activity. This second assumption 
may be valid at high redshift z>2, but may not be at low 
redshift where fueling must be an issue. We therefore com- 
pare their models to our data for z = 2 only and use our 
two bins at z = 1.92 and z = 2.13 to make the compar- 
isons. To convert from £(20) to n we assume an under- 
lying CDM power spectrum with the WMAP/2dF param- 
eters. This results in n(z = 1.92) = 9.35±1;69 hT x Mpc 

h^ 1 Mpc. We also need 



+ 1.58 



and n(z = 2.13) = 11.29 
to convert between the space density assumed by Martini 
& Weinberg (5.27 x 10~ 7 h 3 Mpc" 3 for ft A = 0.7 and 
Sl m = 0.3) and the measured space density of the 2QZ at 
z = 2 (5.1 x 10~ 6 h 3 Mpc~ 3 for the same cosmology). This 
difference increases the estimated lifetimes by a factor of 9.7 
compared to those derived by Martini & Weinberg. We then 
use the Martini & Weinberg fitting function for lifetimes in 
a ACDM Universe (cr 8 = 0.9) to find that t Q = 9.7±{j;| x 10 6 
years (for the z = 1.93 point) and tQ = 2.4+14 x 10 7 years 
(for the z = 2.13 point). Thus the full range of lifetimes at 
z — 2 in this model is £q ~ 4 — 50 Myr. This range is lower 
than, but consistent with the upper limits derived above. 

The above determination of the typical QSO lifetime is 
the total period of activity for a single BH, which may be 
split up into several episodes of activity. The short lifetime 
indicates that there are many generations of QSOs, and that 
a large fraction of galaxies pass through an AGN phase. The 
models used by Martini & Weinberg and others generally as- 
sume that luminosity is perfectly correlated with host mass, 
thus more luminous QSOs would be in more massive DMHs 
and therefore be more strongly clustered. We will investigate 
this below (see Section 6). A scatter in the relation between 
DMH mass and QSO luminosity, would tend to increase the 
effective lifetime, and thus the estimates from the Martini & 
Weinberg models become lower limits to the QSO lifetime. 



5.3.2 Accretion efficiency and the mass of black holes 

There is strong evidence for a correlation between bulge 
velocity dispersion, o~ c and central BH mass (Gebhardt et 
al.2000; Ferrarese & Merritt 2000). This has been extended 
to a correlation between Mbh and Mdh by Ferrarese (2002). 
The exact connection is uncertain, largely due to uncertainty 
in the DMH density profile. Ferrarese suggests three possible 
relations, covering the likely range of allowable assumptions: 



10 8 M q 



0.027 



I 10 12 M o 



(22) 




Figure 24. a) The estimated Mbh based on the relations of Fer- 
rarese (2002) (points connected by solid lines) and Wyithe & Loeb 
(2004) (points connected by dotted lines). We show estimates of 
Mbh based on Eq. 22 (rilled circles) Eq. 23 (filled squares) and 
Eq. 24 (filled triangles) for the Ferrarese (2002) relations and for 
e = £ SIS (open circles), 3.7esiS (open squares) and 25esiS (open 
triangles), b) The derived accretion efficiency, L/Z/Edd; from the 
above Mbh estimates, using the same symbols at in the plot of 
M BH - 
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for an NFW profile (Navarro et al. 1997) and 
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(23) 



(24) 



for a profile based on the weak lensing results of Seljak 
(2002) (henceforth S02). If we assume that these relations do 
not evolve with redshift, then we can directly estimate the 
central BH mass of the DMHs hosting the 2QZ QSOs. These 
BH mass estimates are shown in Fig. 24a (points connected 
by solid lines). We assume h — 0.71 in order to convert from 
/i~ 1 Mq to Mq. As a comparison we also plot estimates of 
Mbh assuming the model of Wyithe & Loeb (2004) in which 
it is the relation between velocity dispersion (or circular ve- 
locity) and Mbh , Mbh — o~ c , which is constant with redshift 
(Shields et al. 2003). This results in a relation between Mdh 
and Mbh of the form 



Mbh = eM D H 



M DH 
10 12 Mr. 



2/3 



A c » m (0) 
187T 2 ft m (z) 



5/6 



(l+ 2 ) 5/2 ,(25) 



where e is a constant and 



18tt + 82(Sl m (z) - 1) - 39(n m («) - l) 2 . 



(26) 



for an isothermal dark matter profile, 



The constant e depends on the density profile of the DMH 
and based on the work of Ferrarese (2002) Wyithe & Loeb 
suggest that for the assumption of a singular isothermal 
sphere e S is ^ 10~ 51 . For a NFW profile e = 3.7esis and 
for an S02 profile e = 25esis- These models with, e = esis, 
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3.7esis and 25esis (which are direct analogues of Eqs. 22, 23 
and 24 for the case of a non-evolving Mbh — c c ) are plotted 
in Fig. 24a (points connected by dotted lines). Examination 
of this plot shows that models in which the Mbh — cr c is 
independent of redshift predict higher mass BHs, and a sig- 
nificant increase in Mbh with redshift for 2QZ QSOs. The 
masses in this case are a factor ~ 50— 100 greater at z = 2.5 
than they are at z = 0.5. In contrast, for the assumption 
that Mbh — Mdh is independent of redshift, there is a much 
weaker trend of increasing Mbh- 

Given the known mean absolute magnitude of each 
redshift interval, we can then calculate the accretion effi- 
ciency, L/LEdd, where L is the bolometric luminosity of 
the QSOs and Z/Edd is the Eddington luminosity [LEdd = 
10 39 1 (M/10 8 M©)W]. To determine the bolometric luminos- 
ity we convert from absolute magnitude in the 6j band using 
the relation derived by McLure & Dunlop (2004) for the B 
band and correcting by bj — B — 0.06 for a mean QSO 
B — V = 0.22 (Cristiani & Vio 1990). The relation is then 

M bj = -2.66 log(L) + 79.42 (27) 

for L in Watts. The resulting accretion efficiencies are shown 
in Fig. 24b. In some cases the mean efficiency of the popula- 
tion is found to be super-Eddington. If the Eddington limit 
is a meaningful constraint on the accretion of matter onto 
super-massive BHs, then the Mbh — Mdh relations described 
by Eqs. 22 and 23 are unlikely to hold at high redshift, 
as they predict that accretion that is significantly super- 
Eddington. For the relation described by Eq. 24, L/Lsdd 
evolves little and is at ~ 0.1 at all redshifts. There is also 
little evidence of evolution for the cases in which Mbh — cr c 
is independent of redshift (connected by dotted lines). The 
values for L/L E dd range between L/L B dd ~ 1 and ~ 0.01 de- 
pending on the value of e assumed. The more realistic values 
of e (3.7esis and 25esis) imply a lower accretion efficiency. 
We note that Wyithe & Loeb (2004) have fit models to the 
QSO clustering results presented by Croom et al. (2001a). 
They suggest that a model where Mbh — o"c is independent 
of redshift is preferred from this data, however, this assumes 
that the accretion efficiency is not a function of redshift. 

An independent estimate of Mbh is available by invok- 
ing the virial theorem in the QSO broad line region and 
using the widths of broad lines as a direct probe of the kine- 
matics. Authors have carried out this analysis on both the 
2QZ (Corbett et al. 2003) and SDSS (McLure & Dunlop 
2004). There are a number of assumptions in these analysis. 
The most crucial of which is the radius-luminosity relation 
for broad line regions (Kaspi et al. 2000). This is generally 
assumed to be independent of redshift, although this has not 
been demonstrated observationally. These works provide a 
relatively independent comparison to the present analysis. 
Corbett et al. (2003) find little evidence of any evolution 
of L/L Ed d in the 2QZ. McLure & Dunlop (2004), also find 
only weak evolution in L/LEdd for the SDSS. Note that both 
of these samples are flux limited so that higher luminosity 
QSOs are at higher redshift, however, it is then still true 
that QSOs with L ~ Lq have little evolution in L/L E dd- 

This implies that the evolution in luminosity of Lq 
QSOs is not caused by a decline in fuelling, but rather, by 
less massive BHs becoming active at lower redshift. It is also 
possible that the observed break in the QSO LF (see Paper 
XII) may be due to the difficulty of accreting with an effi- 



ciency above some limit (presumably close to the Eddington 
limit). However, the shape of the QSO LF is likely driven 
by a combination of accretion rate and Mbh- Any spread in 
accretion rate for a given Mbh would suppress any luminos- 
ity dependence of QSO clustering. We will investigate this 
issue in the next Section. 



6 THE LUMINOSITY DEPENDENCE OF QSO 
CLUSTERING 

In this section we investigate whether there is any evidence 
for QSO clustering being dependent on luminosity. There 
is evidence that low redshift AGN have nuclear luminosities 
that are correlated with host galaxy luminosity (e.g. Schade, 
Boyle & Letawsky 2000), and in particular with the lumi- 
nosity of the bulge/spheroid component of the host. It has 
also been shown that galaxy clustering is a strong function 
of luminosity brighter than L* al (e.g. Norberg et al. 2001). 
Thus bright QSOs, which would be expected to inhabit the 
most massive galaxies, should be clustered more strongly 
that faint QSOs. Croom et al. (2002) investigated this in 
the first data release of the 2QZ (Croom et al. 2001b), and 
found some weak evidence for QSOs with brighter apparent 
magnitudes (approximately equivalent to luminosity relative 
to Lq) being more strongly clustered. A range of physical 
affects could act to cancel any correlation of clustering with 
luminosity. For example, a broad range of accretion efficien- 
cies. 

It is possible to examine the luminosity dependence of 
QSO clustering in a number of ways. Ideally, we would split 
the sample up into a number of redshift and luminosity bins 
and try to separate the luminosity and redshift dependen- 
cies. This is hard simply due to the low number density 
of QSOs, particular in the most luminous intervals. In the 
analysis below we follow Croom et al. (2002) and measure 
the clustering of QSOs as a function of apparent magnitude. 
This has a number of advantages, as it allows us to split the 
QSOs up into only a small number of sub-samples. Apparent 
magnitude is also approximately equivalent to a magnitude 
relative to Lq over the redshift range we are considering, due 
to the strong evolution of the QSO LF. This means that in a 
given apparent magnitude interval, QSOs will have approx- 
imately the same space density at every epoch. 



6.1 QSO clustering as a function of 6j 

We split the 2QZ QSOs into five sub-samples, on the basis 
of their apparent magnitude, bj. These intervals are listed 
in Table 4. To enhance the dynamic range of this analy- 
sis we also include QSOs from the 6dF QSO Redshift Sur- 
vey (6QZ; Paper XII). This data set contains 275 QSOs at 
0.3 < 2 < 2.2 in the magnitude range 16.0 < 6,j < 18.25 se- 
lected from the same photometric data as the 2QZ. It forms 
a bright extension to the 2QZ, in the SGP region only (see 
Paper XII). All the QSOs in the 6QZ form a sixth mag- 
nitude interval. The distribution of QSOs in the z — M bj 
plane is shown in Fig. 25. Even with the large sample pre- 
sented here, the steep bright-end slope of the QSO lumi- 
nosity function means that we can only cover an effective 
dynamic range of ~ 3 mag in apparent magnitude (or a 
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Table 4. 2QZ/6QZ clustering results as a function of apparent magnitude, bj, for a WMAP/2dF cosmology. All fits are over scales 
s = 1 — 25 h^ 1 Mpc. We list the 6j interval and the mean redshift, apparent magnitude and absolute magnitude for each bin together 
with the number of QSOs used. We also give the value of at the mean redshift of each sample derived assuming the polynomial 

evolution model of Paper XII. The best fit values of so (in units of h~ 1 Mpc) and 7 are given with their \ 2 values, number of dof, v 
and probability of acceptance, P(< X 2 )- Lastly wc also list the measured values of for s = 20, 30 and 50 h~ 1 Mpc. We do not fit a 
power law to the brightest magnitude bin (6QZ data) as there are to few QSO-QSO pairs to make a reliable fit, we also don't list f(20) 
for this sample, as there are no pairs found on scales < 20 h Mpc. 
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Figure 25. The rcdshift-absolutc magnitude distribution of 2QZ 
(small points) and 6QZ (triangles) QSOs used in our analysis. 
The solid lines denote the apparent magnitude limits applied to 
the data, while the dashed lines show the redshift range used. A 
WMAP/2dF cosmology is assumed. 



factor of ~ 16 in luminosity). There is also only a rela- 
tively small dynamic range in QSO space density, from a 
mean 4.5 x 10 _6 ft 3 Mpc _3 mag _1 at the faintest magnitudes 
to 9.2 x lCn 7 ft 3 Mpc~ mag . The greatest luminosity de- 
pendence might be expected for the brightest QSOs, as these 
are the rarest sources. This is exactly the point at which the 
rarity of QSOs makes clustering measurements most diffi- 
cult. One solution to this problem is to cross-correlate QSOs 
of a given luminosity with QSOs at all other luminosities. 
This approach will be discussed by Loaring et al. (in prepa- 
ration) . 



The measured fej dependent £(s) are shown in Fig. 26. 
At bright magnitudes (Fig. 26a) the small number and low 
space density of QSOs means that no significant signal is 
detected. At fainter magnitudes the data appear reasonably 
consistent with the best fit power law for the full sample 
(dotted lines). We also fit power laws to each bj interval, 
showing the results as the solid lines in Fig. 26. The val- 
ues are also listed in Table 4. The best fit parameters vary 
considerably, but have large errors. Neither the slopes or am- 
plitudes are particularly well constrained. If instead we fix 
7 = 1.2 as found above, we find values of so that are much 
closer to the mean (dashed lines in Fig. 26). We also note 
that the faintest magnitude interval (Fig. 26f) shows more 
structure on large scales than the other samples. It is possi- 
ble that this is the result of increased incompleteness at the 
faint limit of the sample, even though we have taken care 
to correct for magnitude dependent spectroscopic complete- 
ness, as described in Paper XII. Estimation of £(s) using 
the RA-Dec and RA-Dec-z mixing methods described above 
cause some reduction in this excess at large scales but does 
not completely remove it. This suggests that some, but not 
all, of this excess power could be due to residual incomplete- 
ness affects. Bearing this in mind we have checked whether 
any of our results above are affected by removing QSOs in 
the faintest bin from our sample and confirm that they have 
no significant impact on our conclusions. 

In order to use a robust measure of any luminosity de- 
pendence we calculate £(s) in each of the bj intervals (Table 
4), which is plotted in Fig. 27. We confirm that the esti- 
mates of £ are not significantly changed by using the RA- 
Dec mixing method to measure £(s). We find that there 
is no significant evidence for any dependence of clustering 
amplitude with bj (or equivalently luminosity relative to 
Lq). However, given the relatively large errors found [~ 30 
per cent in £(20)] this result does not rule out models for 
which QSO clustering should be dependent on luminosity. 
As pointed out above, the mean space density of our bright- 
est and faintest samples only differs by a factor of ~ 5. If 
this decrease in space density was solely due to higher mass 
(and therefore rarer) halos acting as hosts then this would 
correspond to a factor of ~ 2 increase in Mdh, but only a 
~ 15 per cent increase in bias (or ~ 30 per cent in cluster- 
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Figure 26. The QSO £(s) from the 2QZ/6QZ (filled points) as a function of apparent fej magnitude in 6 intervals from bright (a) to faint 
(f) magnitudes. The best fit power law is shown in each case (solid line) as is the best fit power law for the full sample for comparison 
(dotted line). We also show the best fit power law when fixing 7 to a value of 1.20 (dashed lines). No power law fit is attempted for the 
6QZ data (a). A WMAP/2dF cosmology is assumed. 
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Figure 27. The dependence of £(s) on bj for three different values 
of s = 20, 30 and 50 h _1 Mpc (open circles, filled circles and open 
squares respectively). We do not plot a point at s = 20, for the 
brightest bin, as there arc no QSO pairs found. 

ing amplitude) which is approximately at the level of our 
measurement errors. This suggests that the increase in sen- 
sitivity provided by cross-correlating different QSO samples 
may provide useful constraints on QSO models (see Loaring 
et al. in preparation). 



7 CONCLUSIONS 

We have performed a detailed analysis of the clustering of 
2QZ QSOs in redshift space as described by the two-point 
correlation function. Here we now discuss our conclusions. 

The QSO two-point correlation function, £(s), averaged 
over the redshift range 0.3 < z < 2.2, shows a slope which 
changes as a function of scale, being flatter on small scales 
and steeper on large scales. A power law is an acceptable 
fit on scales less than 25 ft" 1 Mpc in a WMAP/2dF cos- 
mology, the best fit parameters are so = 5.48l" '4g h" 1 Mpc 
and 7 = 1.20tg;$. We demonstrate that QSO clustering on 
scales < 10 h^ 1 Mpc is strongly affected by non-linear z- 
space distortions, caused by redshift errors, shifts in QSO 
broad emission lines and intrinsic peculiar velocities, which 
all contribute similar amounts to the total velocity disper- 
sion, of {w'i} 1/2 ~ 690 km s" 1 . A power law model which 
has been corrected for both linear and non-linear z-space 
distortions is shown to be a good description of the shape 
of £(s). Here we note that in modelling non-linear z-space 
distortions at high redshift it is important to include an ex- 
tra factor of 1 + z in Eq. 11 relative to the version normally 
used. 

On large scales power law clustering is not an appro- 
priate model and we therefore compare the 2QZ £(s) to a 
model CDM £(s) in a WMAP/2dF cosmology (fi m = 0.27, 
0,a = 0.73, o-g = 0.84) accounting for the affects of non- 
linear clustering on small scales and the affects of z-space 



distortions. This model is well matched to the data after 
allowing for a linear bias of &q = 2.02 ± 0.07 at the mean 
redshift of the sample (z = 1.35). The 2QZ £(s) also agrees 
remarkably well with that measured from the low redshift 
galaxies in the 2dFGRS (Hawkins et al. 2003) , in both shape 
and amplitude. While the match in shape is unsurprising 
given that the physics (at least on large scales) prescribing 
the shape should be identical, the match in amplitude is im- 
pressive. Given that 2dFGRS galaxies are unbiased tracers 
of the mass distribution at low redshift (Verdi et al. 2002), 
it appears that the bias of QSOs exactly cancels out the 
growth of density fluctuations, to give a measured cluster- 
ing equivalent to an unbiased population at low redshift. 
As we find evidence for evolution of QSO clustering in a 
WMAP/2dF cosmology, this agreement must be something 
of a coincidence. Also, in an EdS Universe the 2QZ £(s) is 
a factor ~ 2 below the observed 2dFGRS clustering. 

To further investigate these issues, we determine the 
clustering of 2QZ QSOs as a function of redshift. In a 
WMAP/2dF cosmology we find a significant (at the 98 
per cent level) correlation of clustering amplitude with red- 
shift as measured by the integrated correlation function 
within 20 h' 1 Mpc, £(20). Clustering increases with red- 
shift and we find £(20) = 0.263 ± 0.075 at z = 0.53, and 
£(20) = 0.701 ± 0.174 at z = 2.48. In an EdS cosmology we 
find no evidence for evolution. By assuming an underlying 
WMAP/2dF cosmology we are able to directly determine 
the bias of QSOs, which we find to be a strong function of 
redshift. Even if there were no evolution in the measured 
£(20) with redshift, this would still imply a strongly evolv- 
ing QSO bias. At low redshift, 2QZ QSOs appear largely 
unbiased, with 6q(z = 0.53) = 1.13 ± 0.18, while at high 
redshift we find bq(z = 2.48) = 4.24 ± 0.53. A complication 
is that as the 2QZ is a flux limited sample, we are sam- 
pling more luminous QSOs at high redshift. However, the 
strong evolution of the QSO population means that to good 
approximation we are sampling the QSO population at the 
same space density at each redshift, and at the same point 
relative to the evolving break in the luminosity function, 
Lq. It thus appears that Lq QSOs at low redshift should 
be largely unbiased, and clustered similarly to low redshift 
galaxies. This has indeed been seen by Croom et al. (2004) 
who cross-correlate low redshift (and therefore low luminos- 
ity) 2QZ QSOs with 2dFGRS galaxies and find no difference 
in the clustering properties of the two populations [see also 
Wake et al. (2004)]. 

By using the theoretical relation between Mdh and bias 
derived by Sheth et al. (2001) and others, it is possible for us 
to take the measured bias values for 2QZ QSOs and calculate 
the typical masses of their hosts' DMHs. We find that the 
mass of DMHs hosting 2QZ QSOs is approximately constant 
with redshift, with a mean M DH = (3.0±1.6) x lO 12 /i _1 M . 
The fact that the hosts of 2QZ QSOs have the same mass 
at all redshifts demonstrates that they cannot be cosmologi- 
cally long lived, as DMHs tend to grow and accumulate mass 
over time. Based on the formalism of Lacey & Cole (1993) we 
predict that DMHs hosting QSOs at z ~ 2.5 would typically 
have merged into DMHs of mass ~ 6 x 10 14 ft~ x Mq by the 
present, and therefore exist in rich galaxy clusters (although 
they would not generally be active at low redshift). In con- 
trast, the descendents of lower redshift 2QZ QSOs would 
not have had time to form more massive halos, and should 
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exist in either massive galaxies or groups. By extrapolation 
it is suggested that at z — 0, ~ Lq QSOs should also sit in 
Mdh = (3.0 ± 1.6) x lO 12 /i _1 M halos, which is very close 
to the characteristic mass of the Press-Schechter mass func- 
tion, M*(z = 0) ~ 3.5 x 1O 12 /i _1 M . The predicted growth 
of DMH mass by accretion/merging allows us to place up- 
per limits on the lifetime of the QSO population. Low red- 
shift QSOs cannot be the same population of objects as at 
higher redshift if they have masses which are less then the 
mass of the high redshift sources, after accounting for their 
expected growth over time. Therefore calculating the time 
taken to reach the mean QSO host DMH mass plus twice the 
measured rms gives a ~ 2a limit on the lifetime of QSO ac- 
tivity. We find this limit to be £q < 6 x 10 8 years at z = 2.48, 
but weaker at low redshift (3 x 10 9 years at z = 0.53). We 
note that this limit is not based on the measured number 
density of QSOs compared to a Press-Schechter mass func- 
tion (as many other estimates are), but is only constrained 
by the clustering evolution of QSOs. Various authors have 
provided more detailed models in order to constrain the life- 
time of QSO activity. When applied to our data, the model 
of Martini & Weinberg (2001) suggests that z ~ 2 QSOs 
will have lifetimes £q ~ 4 — 50 x 10 6 years. If there is scatter 
in the relation between Mdh and luminosity, then this is an 
effective lower limit on QSO lifetimes. The e-folding time for 
the evolution of Lq is ~ 2 x 10 9 years (Paper XII), much 
less than the ages determined from the Martini & Weinberg 
model, and significantly less than our clustering evolution 
upper limits at high redshift. 

As a next step we determine the central BH mass, Mbh, 
of 2QZ QSOs based on their estimated Mdh - For this we use 
the relations suggested by Ferrarese (2002) to estimate Mbh 
for different assumptions concerning the density profiles of 
the DMHs, and the evolution of the correlation (Wyithe 
& Loeb 2004). A model in which the correlation between 
Mbh and Mdh is unchanging with redshift predicts that 
BH masses should be slightly increasing with redshift, with 
Alog(MaH) — 1.3 ± 1.1 from the lowest to highest redshift. 
The derived BH masses are in the range 1 - 20 x 1O 7 M 
for NFW profiles, or 0.9 - 20 x 10 8 M Q for S02 profiles. The 
Eddington ratio, L/L^dd, is seen to be approximately con- 
stant as a function of redshift when the Mbh — Mdh relation 
is independent of redshift. This is found to be significantly 
greater than 1 if isothermal DMHs are assumed, and approx- 
imately 1 for the NFW profile, while the S02 profile gives 
L/L-Edd ~ 0.1. A model in which it is the Mbh — <? c relation 
which is invariant with redshift gives a much stronger evolu- 
tion of Mbh as DMHs of a given mass have a higher central 
velocity dispersion when formed at higher redshift. Thus the 
change in Mbh from low to high redshift is more significant 
with Alog(M BH ) ^ 2.1 ± 1.1, and BHs of order ~ 1O 1O M 
being predicted at high redshift. This increase in estimated 
Mbh is greater than (although not significantly) the factor 
~ 30 increase in mean luminosity from our lowest to high- 
est redshift interval. As a result there is a small (factor of 
a few) decline in L/L^dd with increasing redshift, although 
again this is not significant. As the BH masses predicted 
are higher, the accretion efficiencies are lower, in the range 
L/L-Edd ~ 0.01 — 1 depending on DMH profile assumed. 

The above suggests that any model of BH formation 
in which super-massive BHs form at least as efficiently at 
high redshift as they do at low redshift, will tend to have 



L/L^dd constant or decreasing with redshift. This implies 
that it cannot be a reduction in efficiency which is driving 
the fading of the QSO population to low redshift. Instead ac- 
tive BHs at high redshift are more massive that those at low 
redshift, and it is this reduction in the BH mass that causes 
the population of bright QSOs to disappear in the local uni- 
verse. Because super-massive BHs cannot be destroyed, the 
massive BHs active at high redshift must be largely inactive 
at low redshift, otherwise we would find that low redshift 
QSOs would show lower accretion efficiency, and be located 
in more massive DMHs. This argument also implies that at 
any given redshift, the QSO population must be dominated 
by objects which are active for the first time. Hence it is 
likely that each QSO passes through only one bright active 
epoch (possible at the point of BH formation), although at 
low redshift massive BHs may accrete at levels well below 
i/Edd without contributing significantly to the total luminos- 
ity of the population [see also the discussions in Corbett et 
al. (2003) and Croom et al. (2004b)]. 

The above is valid at redshift below z ~ 2.5, which is 
approximately the point at which the space density of lu- 
minous QSOs peaks. Clustering measurements of QSOs at 
z > 2.5 would help us to understand the build up of QSOs at 
this epoch. However the low surface density of z > 2.5 QSOs 
currently makes any accurate clustering measurements dif- 
ficult or impossible. The increasing number of high redshift 
QSOs from the SDSS survey (Fan et al. 2001) may remedy 
this situation. 

Finally, we examine our sample to look for any indi- 
cation of luminosity dependence in the clustering of 2QZ 
QSOs, by measuring £(s) as a function of apparent mag- 
nitude. This shows no indication of any luminosity depen- 
dence that might be expected if more luminous QSOs inhab- 
ited more massive DMHs, but the errors are large enough 
that we would not be able to detect reasonable amounts of 
luminosity dependence. More detailed investigation of this 
problem will be presented by Loaring et al. (in preparation) . 
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